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ABSTRACT 


An  underwater  glider  is  a  cost-effective  underwater  unmanned  vehicle  with  high- 
endurance  for  oceanographic  research  or  naval  applications.  Its  navigation  and 
localization  accuracy  are  important  because  these  accuracies  provide  spatiotemporally 
high  resolution  ocean  data  with  saving  energy  and  time.  The  glider,  however,  is  affected 
by  the  ocean  currents  because  of  its  minimal  velocity,  which  is  due  to  its  buoyancy- 
driven  propulsion  system.  It  also  lacks  of  inexpensive  and  efficient  localization  sensors 
during  its  subsurface  mission.  Therefore,  knowing  its  precise  underwater  position  is  a 
challenging  task. 

This  study  attempts  to  develop  a  novel  correction  method  for  estimating  a  glider’s 
optimal  underwater  trajectory.  In  four  steps,  it  compares  the  corrected  trajectories,  which 
are  developed  using  depth-averaged  and  depth-dependent  correction  methods  using  the 
Regional  Navy  Coastal  Ocean  Model  (NCOM).  The  results  suggest  that  the  depth- 
dependent  correction  method  is  more  accurate.  This  study  for  estimating  a  glider’s 
underwater  trajectory  accurately  would  be  beneficial  to  oceanographic  research  and  naval 
applications,  especially  antisubmarine  warfare  (ASW)  such  as  operating  Intelligence, 
Surveillance,  and  Reconnaissance  (ISR);  operating  littoral  ASW;  providing 
communication  networks;  and  supporting  tactical  oceanography. 
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I.  INTRODUCTION 


This  chapter  deals  with  the  background  of  the  thesis  topic,  the  problem  statement, 
hypothesis  and  explanations,  research  questions,  and  the  benefits  of  this  study.  It  also 
describes  the  scope  and  methodology  and  concludes  with  an  outline  of  the  thesis. 

A.  BACKGROUND 

An  underwater  glider  is  a  small,  cost-effective,  and  multifunctional  unmanned 
underwater  vehicle  (UUV)  designed  for  long-duration  missions  such  as  oceanographic 
research,  naval  operations,  and  especially  antisubmarine  warfare  (ASW).  Since  an 
underwater  glider  is  equipped  with  a  wireless  network  and  requires  less  energy  during  its 
missions,  it  would  be  an  effective  platform  for  future  ASW,  such  as  Intelligence, 
Surveillance,  and  Reconnaissance  (ISR),  and  targeting  and/or  identifying  the  enemy. 

ASW  is  a  branch  of  undersea  warfare  (USW)  conducted  with  the  intention  of 
destroying  and  neutralizing  enemy  submarines  by  using  naval  assets  such  as  surface 
ships,  submarines,  and  aircrafts.  ASW’s  operational  objectives  are  to  deter  the  enemy 
from  using  its  submarines  effectively  and  to  allow  friendly  forces  to  operate  throughout 
the  area  of  operations  (AO).  It  is  crucial  for  naval  forces  to  accomplish  maritime 
superiority  and  to  establish  dominance  in  the  underwater  environment  simultaneously 
(DOD  2013,  Chapter  IV). 

ASW  is  a  vital  naval  operation  because  of  the  continuous  development  of 
submarines  and  underwater  weapons.  ASW,  especially,  has  emerged  as  a  more  important 
issue  for  the  Republic  of  Korea  Navy  (ROKN)  since  the  ROK  Ship  Cheonan  was  sunk  by 
the  North  Korea  Navy’s  torpedo  attack  in  2010.  Since  then,  the  ROKN  has  been  required 
to  improve  ASW  capability.  As  part  of  these  improvements,  the  underwater  glider  is 
expected  to  be  a  suitable  platform  for  future  ASW. 
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B. 


PROBLEM  STATEMENT 


An  underwater  glider  could  be  beneficial  to  future  ASW  in  the  ROKN,  but  it  is 
difficult  to  estimate  its  underwater  trajectory  precisely.  To  conduct  ASW  effectively, 
accurate  localization  of  a  glider  is  important  because  observed  ocean  data  or  the  enemy 
target  information  can  be  more  useful  when  the  location  of  the  data  is  precise.  Localizing 
a  glider’s  underwater  position,  however,  is  difficult  and  complex.  The  glider’s  velocity  is 
limited  compared  to  that  of  other  UUVs  due  to  its  buoyancy-driven  propulsion  system. 
Typically,  the  horizontal  velocity  is  about  0.4  ms'1  and  the  vertical  velocity  is  about  0.2 
ms'1,  which  is  the  same  order  of  magnitude  as  the  ocean  currents  (Smith  et  al.  2010a). 
Thus,  the  glider’s  movement  is  strongly  affected  by  the  ocean  currents.  In  addition,  a 
glider  does  not  have  sensors  for  measuring  its  position  and  the  subsurface  ocean  currents 
in  real  time  (Smith  et  al.  2010b),  and  it  can  only  receive  a  Global  Positioning  System 
(GPS)  fixed  position  when  it  is  located  at  sea  surface.  Consequently,  given  that  some 
margin  of  error  is  common  between  the  dead  reckoning  (DR)  position  (P,ir  in  Figure  1) 
and  the  real  surfacing  position  (Pgps  in  Figure  1),  knowing  its  actual  underwater  location 
is  difficult. 


Figure  1.  Schema  of  Glider  Trajectories. 
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While  improving  the  method  for  estimating  a  glider’s  underwater  trajectory  is 
important,  it  is  a  challenging  task  due  to  several  errors  in  the  existing  methods.  As  of 
now,  estimating  a  glider’s  trajectory  is  done  by  linear  interpolation  between  actual 
surfaced  positions  (T3  in  Figure  1)  or  by  improving  the  accuracy  of  its  DR  trajectories  (Ti 
in  Figure  1).  The  fonner  method,  however,  is  not  realistic  because  the  position  is 
estimated  by  dividing  the  distance  between  the  real  GPS  fixed  position  by  transit  time, 
and  the  glider’s  direction  and  speed  remain  constant.  In  other  words,  the  trajectory  looks 
like  a  straight  line.  The  latter  method  has  two  major  types  of  error:  inherent  error  and 
extraneous  error.  Inherent  error  can  be  caused  by  the  glider’s  dynamic,  mechanical,  and 
measurement  errors,  such  as  heading,  roll  and  pitch,  and  pressure.  On  the  other  hand, 
extraneous  error  can  be  caused  by  environmental  factors,  mainly  the  ocean  currents. 
These  errors  have  an  effect  on  the  estimation  of  DR  trajectory. 

C.  HYPOTHESIS  AND  EXPLANATIONS 

The  accuracy  for  estimating  a  glider’s  underwater  trajectory  can  be  increased  by 
correcting  these  errors  with  two  major  models:  the  glider’s  vehicle  model  and  the  ocean 
circulation  model.  The  former  is  the  solution  to  inherent  error,  and  the  latter  is  the 
solution  to  extraneous  error.  Recently,  Wang  et  al.  (2013)  and  Singh  (2014)  studied 
localization  improvement  for  an  underwater  glider  based  on  its  vehicle  model.  Smith  et  al. 
(20 10a, b)  developed  the  algorithm  for  a  glider’s  path  planning  and  trajectory  design  by 
using  ocean  currents  prediction  model  data.  However,  no  localization  method  using  the 
ocean  currents  prediction  model  data  has  been  studied.  Therefore,  this  study  attempts  to 
develop  a  better  method  for  estimating  a  glider’s  optimal  underwater  trajectory  based  on 
the  U.S.  Navy  Coastal  Ocean  Model  (NCOM)  prediction  data,  especially  the  Regional 
NCOM.  Even  though  we  can  exploit  advanced  sensors  as  studied  by  Woithe  et  al.  (2011) 
and  Somers  (2011)  or  additional  reference  UUVs  as  studied  by  Somers  (2011),  these 
methods  require  more  resources  and,  therefore,  more  financial  backing.  Thus,  this  thesis 
focuses  on  the  single  glider  operating  situation  without  any  advantage  of  additional 
equipment. 
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D.  RESEARCH  QUESTIONS 

The  primary  research  questions  of  this  thesis  are:  “How  can  we  estimate  a  glider’s 
optimal  underwater  trajectory  more  accurately?”  and  “How  does  this  study  impact  future 
ASW  in  the  ROKN?” 

As  a  part  of  the  research  process,  the  following  subsidiary  questions  will  be 
discussed: 

•  Why  is  the  precisely  estimated  underwater  trajectory  of  a  glider 
important? 

•  What  is  the  most  influential  factor  in  estimating  the  underwater  trajectory 
of  a  glider? 

•  Which  method  can  be  used  to  correct  a  glider’s  underwater  trajectory? 

•  How  does  the  method  in  this  thesis  compare  to  the  existing  method? 

E.  BENEFITS  OF  STUDY 

This  study  aims  to  provide  a  better  method  for  estimating  a  glider’s  optimal 
underwater  trajectory,  focusing  on  the  importance  of  the  ocean  currents.  The  precise 
estimation  of  a  glider’s  trajectory  would  be  useful  for  ASW.  During  a  mission,  a  glider 
can  collect  ocean  data,  such  as  Conductivity,  Temperature,  and  Depth  (CTD).  Moreover, 
a  glider’s  navigational  accuracy  would  be  increased  with  well-estimated  trajectories,  and 
this  improved  navigational  accuracy  would  help  it  to  get  more  precise  ocean  data  and  to 
target  information  of  the  enemy.  These  precise  and  valuable  ocean  data  are  critical  for 
calibrating  sonar  to  ensure  that  it  provides  the  most  accurate  underwater  environment. 
ASW  would  be  more  effective  with  more  accurate  infonnation  on  an  enemy’s  target.  As 
a  result,  by  collecting  and  using  these  precise  ocean  data  and  target  infonnation,  the  Navy 
would  be  able  to  improve  ASW  capability,  such  as  operating  ISR,  targeting,  and/or 
attacking  the  enemy’s  underwater  forces. 
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F.  SCOPE  AND  METHODOLOGY 

Since  this  study  is  limited,  the  research  will  not  consider  a  glider’s  dynamic 
model,  mechanics,  and  measurement  errors.  For  this  study,  data  on  the  Slocum  Electric 
glider  from  a  Naval  Postgraduate  School  class  field  experiment  and  the  Regional  NCOM 
predictions  at  Monterey  Bay  from  the  Naval  Oceanographic  Office  (NAVOCEANO)  will 
be  used  to  improve  the  method  for  estimating  a  glider’s  optimal  underwater  trajectory. 
Details  of  the  background  and  the  methods  will  be  described  in  Chapter  II  and  Chapter  III. 

For  this  thesis,  generally,  the  methodology  is  as  follows.  First  of  all,  this  thesis 
provides  background  on  currently  used  underwater  gliders,  their  navigation  scheme  and 
limitations,  and  their  DR  algorithm  and  limitations  by  conducting  a  literature  review.  It 
also  describes  existing  methods  for  estimating  a  underwater  trajectory  of  a  glider  with  a 
discussion  of  their  limitations  and  provides  background  information  on  an  underwater 
glider’s  dynamic  model  and  its  limitations. 

Secondly,  this  thesis  analyzes  data  on  the  Slocum  Electric  glider  and  the  Regional 
NCOM  predictions  and  details  methods  that  will  be  developed  in  this  study.  The  method 
development  process  is  divided  into  two  phases,  and  the  second  phase  is  divided  into  four 
steps.  Phase  One  aims  to  prove  the  importance  of  the  ocean  currents  in  estimating  a 
glider’s  trajectory  by  comparing  the  distance  between  Pgps  and  Pcir  {dx)  to  Pgps  and  PnCom 
( d2 )  in  each  dive  (see  Figure  2). 
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Figure  2.  Schematic  Presentation  of  the  Phase  One. 

Pgps  is  obtained  from  the  GPS  fixed  position  information,  and  Pjr  can  be 
computed  by  a  glider’s  DR  algorithm.  PnCom,  however,  can  be  updated  by  adding  position 
data  of  T  i  and  the  calculated  water  displacement  from  the  Regional  NCOM  prediction 
data  which  can  be  computed  by  numerical  integration,  simple  composite  trapezoidal  rule. 
After  calculating  the  distance  difference,  we  will  determine  how  much  the  DRNCom  (To) 
is  improved  than  Ti  by  comparing  absolute  ( cl2  - d] )  and  relative  {[dx  -d2\/d]  )  distance 

differences  between  the  two  trajectories. 

Phase  Two  attempts  to  develop  a  better  method  for  estimating  a  glider’s  optimal 
underwater  trajectory  by  correcting  Ti  trajectory  with  two  different  methods:  depth- 
averaged  correction  and  depth-dependent  correction.  In  both  methods,  the  glider’s 
velocity  is  relative  to  the  water  velocity  and  its  underwater  position  is  calculated  by  the 
DR  algorithm.  Therefore,  the  only  aspect  that  needs  to  be  calculated  for  estimating  the 
glider’s  underwater  trajectory  is  the  displacement  of  ocean  water.  By  adding  the 
calculated  water  displacement  to  the  glider’s  DR  position  data  and  correcting  errors 
against  Ti  the  corrected  trajectory  will  be  estimated  more  precisely.  To  calculate  the 
glider’s  underwater  position,  it  is  assumed  that  the  T3  is  a  real  estimated  trajectory  in 
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every  dive  because  knowing  a  glider’s  actual  underwater  position  is  difficult  and  is 
limited.  So  the  T3  is  used  as  an  initial  reference  trajectory  to  match  the  depth  and  time- 
consequent  location  of  the  Regional  NCOM  prediction  data  with  other  trajectories  in 
every  correction  method.  The  T3,  however,  is  not  a  realistic  trajectory.  Thus,  to  compare 
precision  of  the  depth-averaged  correction  method  and  the  depth-dependent  correction 
method,  an  optimal  trajectory  will  be  estimated  by  using  iteration  method  from  the 
corrected  trajectory  with  the  depth-dependent  correction  method. 

The  depth-averaged  correction  method  is  similar  to  the  depth-averaged  currents 
calculation  method  (see  Chapter  II).  This  method  considers  the  averaged  ocean  currents 
during  one  dive.  The  water  displacement  at  each  time  ( tn  )  is  calculated  linearly  by 

dividing  the  total  water  displacement  (AX’1  ,A Yw)  by  time  ratio  at  each  time  ( tn/tN  )  (see 
Figure  3). 

The  depth-dependent  correction  method,  however,  is  a  new  approach  to  correct 
error  of  trajectory.  This  method  considers  the  ocean  currents  from  the  Regional  NCOM 
prediction  data,  which  have  four-dimensional  (4-D)  variables,  i.e.,  U ncom(x,  y,  z,t)  : 
longitude  and  latitude,  depth,  and  time,  respectively.  Therefore,  the  more  precise  water 
displacement  can  be  calculated  by  considering  the  depth  and  time-consequent  location 
and  by  using  the  simple  trapezoidal  rule  (see  Figure  4). 
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Figure  3.  Downward  Looking  View  of  Glider’s  Trajectories  with  the  Depth- 

averaged  Correction  Method. 


Figure  4.  Lateral  View  of  the  Eastward  Component  of  Glider’s  Trajectories 
with  the  Depth-dependent  Correction  Method. 


The  corrected  trajectories  are  compared  to  each  other  to  develop  a  better  method 
in  the  following  four  steps,  and  more  details  of  the  methods  in  each  phase  and  step  will 
be  discussed  in  Chapter  III. 
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•  In  Step  One,  Ti  trajectory  is  corrected  by  using  the  depth-averaged 
correction  method.  The  corrected  trajectory  is  referred  to  as  Ci.  In  this  step, 
Ci  cannot  be  compared  to  T3  because  T3  is  not  a  real  trajectory.  Thus,  only 
results  of  corrected  trajectory  will  be  shown. 

•  In  Step  Two,  the  Ti  trajectory  is  corrected  by  using  the  depth-dependent 
correction  method.  The  corrected  trajectory  is  referred  to  as  C2.  In  this  step, 
the  C2  also  cannot  be  compared  to  T3  for  the  same  reason  as  given  in  Step 
One.  Thus,  only  results  of  the  corrected  trajectory  will  be  shown. 

•  In  Step  Three,  an  optimal  trajectory  is  estimated  by  using  the  iteration 
method  with  C2.  This  optimal  trajectory  is  referred  to  as  C3,  which  is 
expected  to  be  closest  to  a  real  trajectory. 

•  In  Step  Four,  two  corrected  trajectories  are  compared  to  an  optimal 
trajectory  (Ci  and  C3,  C2  and  C3)  to  detennine  which  method  is  developed 
well. 

Finally,  this  thesis  discusses  the  results  by  comparing  these  corrected  trajectories. 
Then,  it  details  the  relevance  of  findings  to  future  ASW  in  the  ROKN. 

G.  THESIS  OUTLINE 

The  rest  of  this  thesis  is  divided  into  four  chapters.  Chapter  II  provides 
background  from  a  literature  review.  In  Chapter  III,  a  method  development  process  for 
estimating  a  glider’s  optimal  underwater  trajectory  will  be  detailed,  followed  by  a 
description  of  the  data  set  and  field  experiment  information.  Chapter  IV  analyzes  the  data 
and  discusses  the  results  of  this  study,  and  then  it  describes  the  relevance  of  these 
findings  to  the  navy.  Chapter  V  concludes  the  thesis  and  recommends  future  research 
areas. 
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II.  BACKGROUND 


A  more  cost-effective  and  accurate  localization  method  needs  to  be  developed  for 
the  navigation  systems  of  an  underwater  glider.  A  glider’s  navigation  and  localization 
accuracy  are  important.  However,  since  an  underwater  glider  is  vulnerable  to  the  ocean 
currents  and  lacks  inexpensive  and  efficient  localization  sensors  (Singh  2014),  estimating 
its  precise  underwater  position  is  challenging.  For  these  reasons,  many  researches  and 
studies  have  been  done  to  develop  a  better  method  for  increasing  the  accuracy  of  a 
glider’s  underwater  position  in  recent  decades. 

Although  the  developed  methods  could  get  a  more  accurate  underwater  position 
than  the  traditional  dead  reckoning  (DR)  algorithm,  they  still  have  some  limitations. 
Currently,  three  types  of  underwater  glider  are  widely  used:  the  Slocum,  the  Spray,  and 
the  Seaglider.  They  have  roughly  similar  appearance  features,  but  their  navigation 
algorithms  and  localization  methods  are  slightly  different.  Thus,  their  navigation 
algorithms  and  localization  methods  must  be  considered  when  generating  a  novel  method 
for  precise  localization.  Most  of  the  previous  studies  developed  methods  by  using  the 
glider’s  vehicle  model,  and  some  studies  measured  the  ocean  currents  with  extra  onboard 
sensors  or  localized  the  glider’s  position  with  additional  equipment.  The  vehicle  model- 
based  methods,  however,  cannot  fully  consider  the  ocean  current  effects  because  they 
focus  on  improving  the  glider’s  attitude,  state,  and  control.  The  methods  using  extra 
onboard  sensors  can  measure  the  ocean  currents  in  near  real-time  and  compensate  for 
error  in  improving  a  glider’s  navigation  accuracy,  but  these  sensors  make  a  glider  heavier 
and  require  more  energy.  Consequently,  the  glider’s  operation  duration  is  shorter  than  it 
was  before.  The  methods  using  multiple  platforms  like  reference  UUVs  and  acoustic- 
based  network  sensors  (Pauli  et  al.  2014)  can  localize  the  glider  in  near  real-time,  but 
they  are  expensive,  and  some  equipment  may  be  subjected  to  the  ocean  environment;  for 
example,  the  ocean  currents  and  bottom  topography  (Pauli  et  al.  2014;  Wang  et  al.  2013). 
Therefore,  this  study  attempts  to  develop  a  more  cost-effective  and  reliable  method  based 
on  the  ocean  model. 
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The  rest  of  this  chapter  explains  the  motivations  for  developing  a  novel  method 
for  estimating  a  glider’s  underwater  trajectory  by  reviewing  the  glider’s  6  Degrees  of 
Freedom  (6  DOF)  non-linear  dynamic  equations  of  motion  and  the  existing  model-based 
methods.  The  following  section  begins  by  presenting  general  descriptions  of  the  three 
commonly-used  underwater  gliders,  including  their  common  navigation  scheme  and  DR 
algorithm.  Then,  each  glider’s  navigation  algorithm  and  its  limitations  are  detailed  to 
provide  background  knowledge  for  developing  a  better  method  for  estimating  a  glider’s 
underwater  position. 

A.  UNDERWATER  GLIDERS 

At  the  present  time,  three  major  types  of  underwater  glider  are  commonly  used: 
the  Slocum  (see  Figure  5a)  that  was  produced  by  the  Teledyne  Webb  Research,  the  Spray 
(see  Figure  5b)  that  was  co-developed  by  Scripps  Institution  of  Oceanography  (SIO)  and 
Woods  Hole  Oceanographic  Institution  (WHOI),  and  the  Seaglider  (see  Figure  5c)  that 
was  developed  by  the  University  of  Washington.  Although  these  gliders  have  roughly 
similar  appearance  and  operational  characteristics,  their  navigational  characteristics  are 
slightly  different.  Such  differences  must  be  studied  and  considered  to  improve  the 
method  for  estimating  a  glider’s  optimal  underwater  trajectory. 


Figure  5.  Three  Types  of  Underwater  Glider:  (a)  the  Slocum  (from 
Autonomous  Undersea  Vehicle  Application  Center  [AUVAC]  2014, 
http://auvac.org/uploads/configuration/Slocuml.jpg);  (b)  the  Spray 
(from  AUVAC  2014,  http://auvac.org/uploads/configuration/spray.jpg); 

(c)  the  Seaglider  (from  University  of  Washington  2014, 
http://www.washington.edu/news/files/20 1 3/05/glider-500x33 1 1  .jpg). 
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1. 


General  Descriptions 


Commonly-used  underwater  gliders  are  small,  cost-effective,  and  multifunctional 
UUV  designed  for  high-endurance  missions.  They  are  approximately  two  meters  long, 
weigh  about  52  kilograms,  and  have  a  hull  diameter  of  about  20-30  centimeters,  so  they 
are  small  enough  to  be  handled  by  one  to  three  persons  (Davis  et  al.  2002).  Due  to  the 
lack  of  propulsion  system,  they  glide  through  the  ocean  by  controlling  their  buoyancy 
using  interior  pumps  and  tanks  and  by  controlling  their  attitude  using  two  wings  and  a  tail 
that  are  attached  to  their  bodies  (Graver  2005).  Depending  on  their  type  and  the  purposed 
of  their  mission,  gliders  are  designed  to  operate  at  the  maximum  depth  of  200-1,500 
meters.  The  Spray  and  the  Seaglider  are  designed  for  deep,  open-ocean  operations,  so 
they  can  dive  at  the  maximum  depth  of  1,000-1,500  meters.  The  Slocum  Electric,  on  the 
other  hand,  is  optimized  for  shallow-water  and  coastal  operations,  and  its  maximum 
depth  is  200  meters.  The  underwater  gliders  move  forward  at  a  horizontal  average 
velocity  of  about  0.4  ms'1  with  a  saw-tooth  vertical  profile  (see  Figure  6).  This  forward 
movement  is  produced  by  changing  the  gliders’  water  displacement  with  a  buoyancy 
engine  (Schofield  et  al.  2007)  and  by  converting  their  buoyancy  into  horizontal 
movement  by  a  lift  force,  which  is  created  from  the  gliders’  wings  (Bender  et  al.  2008). 


Surfacing  Interval 

10— 15min 

Yo-yo  Profile 

3— 12h 

Figure  6.  The  Slocum’s  Saw-tooth  Profile  (from  Flemandez  et  al.  2014). 
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A  glider  is  more  cost-effective  than  other  UUVs  and  ships  (Graver  2005).  The 
glider’s  buoyancy-driven  propulsion  system  allows  it  to  have  minimal  velocity  and  thus 
expends  less  power  than  other  kinds  of  UUVs.  Therefore,  it  can  observe  and  collect  vast 
ocean  data  autonomously  while  it  is  conducting  long-range  and  long-duration  missions. 
In  addition,  a  glider  can  receive  commands  remotely  and  transmit  collected  data  in  near 
real-time  through  a  wireless  telecommunication  systems,  such  as  Radio  Frequency  (RF) 
modem.  Iridium  satellite,  Advanced  Research  and  Global  Observation  Satellite  (ARGOS), 
and  tele-sonar  modem  (Webb  et  al.  2001)  via  antennas  that  are  attached  their  body.  For 
these  reasons,  the  underwater  glider  can  be  used  in  various  missions  like  ocean  sampling 
and  monitoring  (Schofield  et  al.  2007),  as  well  as  naval  applications  such  as  ASW, 
maritime  surveillance  and  reconnaissance. 

a.  Common  Navigation  Scheme 

Although  each  glider  has  unique  characteristics,  the  navigation  scheme  of  the 
three  widely-used  gliders  is  roughly  similar.  Their  navigation  primarily  relies  on  a  GPS 
signal,  magnetic  compass,  altimeter,  and  subsurface  DR.  When  a  glider  locates  at  surface, 
its  actual  position  is  fairly  precise  because  it  can  receive  GPS  fixed  information.  During 
the  subsurface  mission,  however,  its  position  information  is  not  as  accurate  as  the  former 
situation  due  to  the  lack  of  GPS  signals.  Once  the  glider  is  diving  into  the  ocean,  it  is 
gliding  vertically  with  several  upward  and  downward  dives  (see  Figure  6),  until  it 
resurfaces  (Merckelbach  et  al.  2008).  During  the  dives,  the  glider  measures  parameters 
such  as  heading,  pitch  and  roll,  and  depth  change  rate  from  the  onboard  sensors,  i.e., 
attitude  and  pressure  sensors  (Singh  2014).  Then,  the  glider  estimates  its  underwater 
position  by  computing  the  parameters  with  the  DR  algorithm;  the  details  of  the  DR 
algorithm  will  be  described  in  the  next  section.  After  several  dives,  a  glider’s  DR  position 
(see  position  la  in  Figure  7)  is  different  from  the  actual  measured  surfacing  position  (see 
position  1  in  Figure  7).  This  difference  occurs  due  to  the  ocean  currents  and  is  commonly 
referred  to  as  depth-averaged  currents.  The  magnitude  of  depth-averaged  currents  can  be 
calculated  from  the  distance  between  the  DR  and  actual  surfacing  positions,  and  the 
distance  is  divided  by  the  subsurface  transit  time  (Merckelbach  et  al.  2008). 
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A  glider  uses  the  depth-averaged  currents  by  means  of  correcting  the  next 
waypoint  to  compensate  for  these  errors,  and  this  increases  its  navigation  accuracy. 
Merckelbach  et  al.  (2008)  describe  the  Slocum’s  currents  correction  algorithm  (see 
Figure  7),  which  is  similar  to  the  other  gliders,  and  the  details  of  each  stage  are  described 
as  follows. 


underwater  at  surface 

to  t\  ■  t?  <3 

_Z _ Z _ Z _ Z_ 

pre-dive  GPS  resurface  post-dive  GPS  pre-dive  GPS 


Figure  7.  Schema  of  the  Slocum’s  Currents  Correction  Algorithm  (from 

Merckelbach  et  al.  2008). 
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•  At  stage  0,  the  Slocum  prepares  to  dive  underwater  and  obtains  GPS  fixed 
information. 

•  At  stage  1,  the  Slocum  resurfaces  and  prepares  to  receive  GPS  signals  and 
commands  for  the  next  mission  via  telemetry  communications.  Note  that 
after  it  has  resurfaced,  the  Slocum  is  affected  by  surface  currents  until  it 
receives  GPS  signals.  So,  its  actual  position  cannot  be  calculated 
immediately  but  can  be  estimated.  First,  the  surface  current  velocity  can 
be  estimated  from  the  distance  between  position  2  and  3,  divided  by  time 
between  t2  and  t3.  And  then,  position  1  can  be  estimated  from  position  2  by 
using  the  surface  current  velocity  with  the  elapsed  time  between  fi  and  t2. 

•  At  stage  la,  the  Slocum’s  resurfacing  location  is  estimated  by  DR 
algorithm. 

•  At  stage  2,  the  Slocum  obtains  its  first  GPS  fixed  infonnation  after 
resurfacing. 

•  Finally,  the  Slocum  prepares  to  dive  again  and  obtains  GPS  fixed 
information  prior  to  diving. 

b.  Dead  Reckoning  (DR)  Algorithm 

An  underwater  glider  primarily  uses  the  DR  algorithm  to  estimate  its  underwater 
position  and  to  increase  its  navigational  accuracy.  A  glider’s  new  position  is  estimated  by 
the  DR  algorithm  at  about  every  four  seconds  of  the  control  cycle  (Woithe  et  al.  2011) 
using  the  internal  navigation  system.  The  estimated  position  is  calculated  from  the 
glider’s  horizontal  speed  ( Ug  )  and  heading  ( h  )  (Merckelbach  et  al.  2008),  and  the 

horizontal  speed  is  derived  from  the  pitch  angle  ( 0 )  and  depth  change  rate  ( dg )  that  are 

measured  by  the  onboard  attitude  and  pressure  sensors  (Woithe  et  al.  2011).  The  attitude 
sensor  measures  the  glider’s  pitch  angle,  roll,  and  heading,  and  the  pressure  sensor 
measures  their  depth.  Woithe  et  al.  (2011)  present  a  simple  calculation  method  to 
estimate  a  glider’s  underwater  position  and  explain  related  parameters  with  the  following 
equations. 
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The  glider’s  horizontal  speed  (U  )  is  calculated  from  the  depth  change  rate  ( dg ) 

and  present  pitch  angle  ( 6)  during  one  control  cycle  (At)  (Woithe  et  al.  2011),  as  shown 
in  Equation  (1). 


(1) 


Then,  the  glider’s  eastward  ( u  )  and  northward  (  v  )  horizontal  velocity  are 


computed  using  the  glider’s  horizontal  speed,  and  the  measured  present  heading  ( h  ),  as 
shown  in  Equation  (2)  (Merckelbach  et  al.  2008). 


u  =  Ug  sin  h 
v  =  Ug  cos  h 


(2) 


These  components  are  integrated  by  the  real-time  difference  (At)  since  the  last 
control  cycle,  then  converted  to  meters  ( Ax  ,  Ay  )  as  in  Equation  (3). 

Ax-uAt 

(3) 

Ay  —  v  At 

Equation  (4)  represents  the  updated  DR  position  ( x ,  y  )  by  adding  the  newly 
calculated  location  (  Ax ,  Ay  )  with  the  previous  one  (  x ,  y  ). 


x  =  x  +  Ax 

y  =  y  +  4f 


(4) 


Finally,  the  horizontal  distance  ( d  )  during  the  present  control  cycle  time  (At)  is 
calculated  as  shown  Equation  (5)  (Woithe  et  al.  2011). 


Even  the  glider’s  DR  algorithm  is  a  more  accurate  method  for  estimating  the 
underwater  position  than  a  simple  linear  interpolation,  but  errors  between  real  and  DR 
position  exist.  The  main  error  can  be  accounted  for  by  the  ocean  currents.  Some  errors 
are  made  in  the  DR  algorithm  itself  or  are  derived  from  the  simplification  of  glider 
dynamics.  Moreover,  the  others  are  parameter  errors  measured  by  onboard  sensors. 


17 


2. 


The  Slocum 


The  Slocum  is  a  torpedo-shaped  underwater  glider,  which  consists  of  two 
different  models:  the  Slocum  Electric  (see  Figure  5a)  and  the  Slocum  Thennal  (see 
Figure  8).  The  Slocum  Electric  is  propelled  by  a  battery-powered  buoyancy  engine  and  is 
optimized  for  shallow-water  coastal  operation  (Davis  et  al.  2002)  with  various  depth 
ranges  of  4-200  meters,  so  it  is  equipped  with  a  navigation  system  that  can  quickly 
change  the  horizontal  and  vertical  movement  (Webb  et  al.  2001).  It  can  steer  in  a 
horizontal  direction  by  using  a  tail  fin  rudder  (Scholfield  et  al.  2007)  with  a  turning  circle 
of  within  about  seven  meters  (Davis  et  al.  2002).  Also  it  can  rapidly  change  its  buoyancy 
by  pushing  water  in  and  out  through  a  nozzle  that  is  located  in  the  nose  by  using  a  large- 
volume,  single-stroke,  piston  pump  (Davis  et  al.  2002).  The  Slocum  G2,  another  type  of 
the  Slocum  Electric,  was  developed  recently.  It  can  optionally  be  equipped  with  various 
onboard  sensors  which  are  already  developed,  and  it  is  capable  of  operating  with  an 
additional  depth  of  350  meters,  and  1,000  meters  by  changing  some  components  like  the 
buoyancy  engine  (Teledyne  Webb  Research  2014).  As  a  result,  the  Slocum  Electric  has  a 
multi-depth  and  multifunctional  capability.  It,  on  the  other  hand,  would  have  a 
disadvantage  of  much  poorer  applicability  than  other  gliders  when  it  has  the  newly 
developed  sensors  installed.  The  Slocum  Thermal  is  optimized  for  deep-water  operation 
at  a  maximum  depth  of  1,200  meters  and  is  propelled  by  heat  from  the  ocean’s  thermal 
stratification  (Rudnick  et  al.  2004). 
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Figure  8.  The  Slocum  Thermal  (from  Teledyne  Webb  Research  2014, 
http://www.webbresearch.com/thermal.aspx). 

The  Slocum  navigates  using  a  simple  DR  algorithm,  but  this  method  has  some 
limitations.  Both  the  Slocum  Electric  and  Thermal  navigate  between  each  GPS  fixed 
position  by  calculating  their  heading  and  glide  angle,  and  those  parameters  are 
maintained  until  they  are  recalculated  (Bender  et  al.  2008).  During  the  subsurface  mission, 
the  Slocum  calculates  its  buoyancy,  pitch  angle,  and  vertical  velocity  to  estimate  its 
underwater  position  using  the  DR  algorithm  (Webb  et  al.  2001).  The  Slocum,  however, 
does  not  estimate  the  ocean  currents  for  calculating  the  underwater  position  (Bender  et  al. 
2008),  so  its  underwater  trajectory  is  less  precise. 
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Figure  9.  The  Slocum  Electric  Interior  (from  AUVAC  2014,  http://auvac.org/ 
uploads/configuration/Gliderstructure.jpg). 

3.  The  Spray 

Like  the  Slocum,  the  Spray  is  propelled  by  a  battery-powered  buoyancy  engine, 
but  its  body  is  a  nearly  streamlined  torpedo-shaped,  so  the  Spray’s  water  resistance  is 
50%  less  than  the  Slocum  (Sherman  et  al.  2001).  Its  minimal  water  resistance  allows  the 
Spray  to  move  quickly,  using  less  energy  to  overcome  the  strong  water  currents. 
Additionally,  the  Spray  is  less  affected  by  the  strong  ocean  currents  relative  to  the  other 
gliders  because  it  can  control  its  roll  precisely,  thereby  its  attitude  can  be  fixed  well.  It  is 
optimized  for  deep-water  operation  (Rudnick  et  al.  2004)  at  a  maximum  depth  of  1,500 
meters.  The  Spray  controls  its  attitude  by  using  its  internal  batteries:  the  roll  is  controlled 
by  moving  them  to  the  left  or  right,  and  the  pitch  is  controlled  by  moving  them  to  forward 
or  backward  (Sherman  et  al.  2001).  This  control  method,  however,  could  be  inefficient 
because  the  Spray  spends  too  much  time  doing  so  to  change  its  path  when  operating  in  a 
coastal  area  containing  many  obstacles  and  where  the  depth  change  is  complex. 

To  perfonn  the  mission,  the  Spray  rolls  90  degrees  at  the  surface  to  make  one 
wing  vertical  (see  Figure  10)  for  obtaining  GPS  fixes  or  communicating  through  Iridium 
satellite  (Rudnick  et  al.  2004).  After  that,  it  uses  a  simple  navigation  algorithm  to 
compute  the  heading  and  glide  angle,  which  are  required  to  reach  a  desired  location  from 
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its  present  position  with  the  assumption  of  a  constant  pitch  angle,  heading,  and  angle  of 
attack  during  its  subsurface  mission  (Sherman  et  al.  2001). 


Figure  10.  Schema  of  the  Spray’s  Navigation  Feature  (from  SIO  2014, 
http  ://spray.ucsd  .edu/pub/rel/info/ spray_description.php) . 

Like  the  Slocum,  the  Spray  uses  the  DR  algorithm  for  estimating  its  underwater 
trajectory,  but  this  method  is  also  less  precise  because  it  does  not  estimate  the  ocean 
currents.  Therefore,  it  is  difficult  to  navigate  by  following  the  preplanned  path.  The  Spray, 
however,  can  compensate  for  some  errors  on  its  path  by  calculating  the  depth-averaged 
ocean  currents  information. 

4.  The  Seaglider 

The  Seaglider  is  designed  for  deep-water  missions  like  the  Spray.  Its  fairing  is 
enclosed  in  low-drag-shaped  compressible  fiberglass  (Eriksen  et  al.  2001),  which  allows 
the  glider  to  prevent  loss  of  buoyancy  resulting  from  its  depth  change.  Thus,  it  can  dive 
into  deep-ocean  by  using  less  energy  than  other  gliders  which  have  an  incompressible 
body  (Davis  et  al.  2002).  On  the  other  hand,  it  is  difficult  to  load  additional  equipment 
near  the  glider’s  body  because  the  contraction  and  expansion  of  the  body  are  repeated 
depending  on  the  pressure  during  the  subsurface  mission.  Moreover,  the  Seaglider  can 
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only  accommodate  additional  sensors  near  its  wings  due  to  its  structural  features,  so  it  has 
less  versatility  than  other  gliders  (Eriksen  et  al.  2001). 

Unlike  the  Slocum  and  the  Spray,  the  Seaglider  has  a  unique  ability  to  estimate 
the  ocean  currents  by  using  a  Kalman  Filter  (KF)  algorithm  to  increase  its  navigational 
accuracy  (Eriksen  et  al.  2001).  The  Seaglider  automatically  controls  its  direction  to 
compensate  for  the  navigational  errors  between  the  preplanned  path  and  actual  surfaced 
position  based  on  the  DR  position  and  the  indirectly  measured  speed  of  ocean  currents  by 
using  the  KF.  By  using  the  filter,  the  Seaglider  can  approach  the  prescribed  position  more 
precisely  and  expend  less  energy  as  well  (Eriksen  et  al.  2001). 

B.  MODEL-BASED  ESTIMATION  OF  A  GLIDER’S  UNDERWATER 
TRAJECTORY 

Underwater  glider’s  navigation  and  localization  accuracy  are  important  issues 
because  they  allow  a  glider  to  get  spatiotemporally  higher  resolution  data  and  to  save 
energy  consequent  to  its  extended  mission  time.  Precise  localization  of  a  glider  is 
essential  to  obtain  its  navigational  accuracy,  but  it  is  a  challenging  mission.  This  is 
because  a  glider  can  only  receive  GPS  fixed  position  infonnation  when  it  is  surfacing, 
and  GPS  signals  cannot  penetrate  into  the  ocean,  so  a  glider  is  highly  dependent  on  its 
navigation  algorithm  to  estimate  its  underwater  position.  In  addition,  due  to  its  minimal 
velocity  relative  to  the  other  UUVs,  the  attitude  of  a  glider  is  strongly  affected  by  the 
ocean  currents,  and  thus  localization  of  a  glider  is  a  complex  task. 

Most  of  the  underwater  gliders  primarily  use  the  traditional  DR  algorithm  to 
increase  their  navigational  accuracy.  This  method,  however,  still  has  limitations.  First  of 
all,  a  glider  uses  the  depth-averaged  currents  to  compensate  for  error  and  to  set  the  next 
waypoint  (Merckelbach  et  al.  2008),  but  this  method  is  relatively  inaccurate.  Since  the 
depth-averaged  currents  are  estimated  from  the  previous  dive  sequences  in  the  previous 
mission  area,  those  data  are  not  coincident  with  the  areas  of  the  future  mission.  Therefore, 
the  corrected  waypoint  by  the  depth-averaged  currents  would  not  be  able  to  account  for 
the  real-time  ocean  currents  fully.  Second,  the  glider  measures  parameters  by  onboard 
sensors  to  use  them  for  the  DR  algorithm,  but  measurement  errors  still  exist.  These  errors 
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have  a  significant  impact  on  the  magnitude  of  parameters,  which  are  determined  by  the 
glider’s  dynamic  model  and  consequent  localization.  For  these  reasons,  two  types  of 
methods  can  be  considered  to  compensate  for  these  errors.  The  first  method  considers  a 
gilder’s  vehicle  model,  which  determines  the  glider’s  position  and  orientation  in  6  DOF. 
The  six  different  motion  components  of  the  glider  are  measured  by  onboard  sensors  and 
are  calculated  by  its  dynamic  equation  of  motion.  Thus,  the  glider’s  more  accurate 
position  can  be  estimated  by  using  the  well-developed  vehicle  model  and  by  reducing  the 
measurement  error.  The  second  method  considers  the  environmental  model,  especially 
the  ocean  currents,  to  compensate  for  errors  from  external  forces  and  moments  acting  on 
the  glider  (Fossen  1994).  This  model  considers  the  environmental  disturbances — such  as 
ocean  currents,  wind,  and  waves — that  affect  the  glider’s  position  and  motion.  These 
disturbances  can  be  approximated  by  applying  the  principle  of  superposition  and  be 
added  to  the  dynamic  equation  of  motion. 

Some  of  researches  and  studies  (Graver  et  al.  2003;  Graver  2005;  Bender  et  al. 
2008)  discussed  the  dynamic  equations  of  motion  and  parameter  identification  of 
underwater  gliders  in  past  years  (Smith  et  al.  2010b).  This  section  describes  the  6  DOF 
non-linear  dynamic  equations  of  motion  and  discusses  previous  work  for  each  model- 
based  method  to  provide  the  background  knowledge  for  developing  a  novel  method  for 
estimating  a  glider’s  optimal  underwater  trajectory. 

1.  6  DOF  Non-linear  Dynamic  Equations  of  Motion 

A  glider’s  dynamic  model  represents  its  motion  by  a  system  of  equations 
corresponding  to  the  6  DOF,  i.e.,  surge  ( u  ),  sway  (v),  heave  (w),  roll  ( p  ),  pitch  ( q  ), 
yaw  (r),  and  a  glider’s  position  and  orientation  are  detennined  by  the  6  DOF  (see  Table 
1).  The  first  three  DOF  components  represent  the  glider’s  position  and  transitional 
motion,  and  the  last  three  DOF  components  indicate  their  orientation  and  rotational 
motion  (Fossen  1994). 
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Table  1.  Notation  of  6  DOF  for  the  Underwater  Glider  (after  Fossen  1994). 


DOF 

Forces  and 
Moments 

(Body-fixed) 

Linear  and 
Angular 
Velocities 
(Body-fixed) 

Positions  and 
Euler  Angles 

(Earth-fixed) 

1 

Motions  in  the  x-direction  (surge) 

X 

u 

X 

2 

Motions  in  the  y-direction  (sway) 

Y 

V 

y 

3 

Motions  in  the  z-direction  (heave) 

Z 

w 

z 

4 

Rotation  about  the  x-axis  (roll) 

K 

P 

</> 

5 

Rotation  about  the  y-axis  (pitch) 

M 

q 

e 

6 

Rotation  about  the  z-axis  (yaw) 

N 

r 

¥ 

Figure  1 1 .  Earth-fixed  and  Body-fixed  Coordinate  Frames. 

To  develop  the  equations  of  motion,  some  following  assumptions  are  required. 
First,  the  underwater  glider  behaves  as  a  rigid  body.  Second,  the  Earth’s  rotation  can  be 
neglected  because  the  motion  of  the  Earth  rarely  has  an  impact  on  low-speed  gliders. 
Thus,  the  Earth-fixed  frame  can  be  regarded  as  an  inertial  frame.  With  these  assumptions, 
the  glider’s  motion  is  described  in  two  coordinate  frames:  the  Earth-fixed  frame  and  the 
body- fixed  frame  (see  Figure  11).  The  origin  of  the  Earth- fixed  frame  is  the  center  of  the 
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Earth  (O),  and  the  position  and  Euler  angles  of  the  glider  are  expressed  relative  to  this 
reference  frame.  The  origin  of  the  body-fixed  frame,  on  the  other  hand,  is  the  center  of 
the  glider  (CG),  and  the  linear  and  angular  velocities  are  expressed  in  this  frame,  and  the 
motion  of  the  body-fixed  frame  is  represented  relative  to  the  Earth-fixed  frame  (Fossen 
1994). 

Fossen  (1994)  describes  the  general  motion  of  the  glider  in  6  DOF  as  following 
vectors. 

il  =  [x,  y,  z,  (j),  0,  i//f 
v  =  [u,  v,  w,  p,  q,  rf 
t  =  [X,  Y,Z,K,M,N]t 

Here,  q  indicates  the  position  and  Euler  angles  vector  in  the  Earth-fixed  frame,  v 
means  the  linear  and  angular  velocity  vector  in  the  body-fixed  frame,  and  T  denotes  the 
forces  and  moments  acting  on  the  vehicle  in  the  body-fixed  frame.  With  these  6  DOF,  the 
non-linear  dynamic  equation  of  motion  can  be  obtained  by  considering  the  relationship 
between  two  coordinate  frames  (see  Equation  (6)). 

dx 

M  —  +  C(v)v  +  D(x)x  +  g(fi)  =  x  (6) 

dt 

where 

M  :  inertia  matrix  (including  added  mass) 

C :  matrix  of  Coriolis  and  centripetal  terms  (including  added  mass) 

D :  damping  matrix 

g  :  vector  of  gravitational  forces  and  moments 

t  :  vector  of  control  inputs 

After  that,  the  rigid-body  equation  of  motion  is  derived  by  applying  the 
Newtonian  and  Lagrangian  mechanics,  and  the  equation  of  translational  (see  Equation  (7)) 
and  rotational  motion  (see  Equation  (8))  is  obtained  by  applying  Euler ’s  First  and  Second 
Axioms. 

m[u  - vr  +  wq-xG(q 2  +  r2)  +  yG(pq -r)  +  zG(pr  +  q)]  =X 

m[v-wp  +  ur-yG(r2  +  p2)  +  zG(qr- p)  +  xG(qp  +  r)\  =Y  (7) 

m[w -uq  +  vp  —  zG{p2  +  q2)  +  xG(rp  -  q)  +  yG  ( rq  +  p)\  =  Z 
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(8) 


K  'P  +  (f  - Iy)qr - (r  +  pq)I„  +  (r2  - q2)Iyz  +  (pr - q)fy 

+  m[yG  ( w  —  uq  +  vp )  -  zG  (v  -  wp  +  ur)]  =  K 
!yq  +  (4  -  f:  )'V  ~(P  +  qr)Ixy  +  (p2  -  r2)Izx  +  (qp  - r)Iyz 

+  m[zG  (ii  -  vr  +  wq)  -  xG(w -  uq  +  vp)}  -M 
Ij  +  (ly  -  Ix)pq  —  (q  +  rp)Iyz  +  (q2  -  p2)Ixy  +  (rq  -  p)Izx 

+  m[xG  (v  —  wp  +  ur)  -  vG  (;'/  -  vr  +  wq)\  =  N 

where 

xG ,  yG ,  zG  :  center  of  gravity 

I x ,  I y ,  Iz  :  the  moments  of  inertia  about  the  x,  y,  z  axes 

/  ,  /A_ ,  /  :  the  products  of  inertia 

The  Equation  (7)  and  Equation  (8)  can  be  re-expressed  as  the  vectorial 
representation  (see  Equation  (9))  by  considering  the  glider’s  hydrodynamics  and 
moments  with  Radiation-induced  forces  and  Froude-Kriloff  and  Diffraction  forces 
(Fossen  1994). 

M  —  +  C(v)v  +  .D(v)v  +  g-('i|)  =  te  +t  (9) 

dt 

where 

M  =  M  RB  +  M  A 
C(y)  =  Cm  (\)  +  CA(\) 

T  —  T  +  T  +  T 

E  current  wind  wave 

rE  :  the  environmental  forces  and  moments  acting  on  the  vehicle 

t  :  the  propulsion  forces  and  moments 

The  left-hand  side  of  Equation  (9)  expresses  the  glider’s  vehicle  model,  and  the 
right-hand  side  of  the  equation  describes  the  environmental  effects  through  te  .  More 
details  of  each  equation’s  derivation  and  an  explanation  of  terms  are  described  in  Fossen 
(1994). 
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2. 


Previous  Works 


This  section  describes  previous  works  in  developing  methods  for  estimating  the 
underwater  trajectory  of  a  glider  based  on  each  model. 

a.  Vehicle  Model-Based  Methods 

As  previously  discussed,  an  underwater  glider  estimates  its  subsurface  position 
with  the  DR  algorithm  by  measuring  pitch  angle,  heading,  and  depth  change  rate  from  its 
onboard  sensors.  During  the  long-duration  missions,  measurement  errors  of  these 
parameters  can  cause  inaccuracy  in  estimating  a  glider’s  underwater  DR  position.  Since 
these  measured  parameters  are  derived  from  the  glider’s  vehicle  model,  more  accurate 
estimation  can  be  obtained  by  developing  the  model  (Wang  et  al.  2013). 

Most  of  the  studies  consulted  for  this  thesis  developed  methods  for  localization  of 
a  glider  by  using  its  vehicle  model,  but  this  model  has  some  errors.  Woithe  et  al.  (2011) 
developed  a  method  to  improve  the  accuracy  of  the  traditional  DR  algorithm  by  using  a 
Doppler  Velocity  Log  (DVL).  DVL  uses  four  or  more  beams  of  acoustic  wave  to  track 
the  seabed,  and  the  glider’s  relative  motion  (u,v,w)  is  detennined  by  measuring  the 
Doppler-shifted  returns  from  the  seabed  (Pauli  et  al.  2014).  Then,  the  glider’s  position  is 
estimated  by  replacing  the  calculated  velocity  with  traditionally  calculated  speeds  for  the 
DR  process  (Woithe  et  al.  2011).  This  method  can  estimate  the  glider’s  underwater 
position  more  precisely,  but  it  has  some  drawbacks.  Since  the  extra  onboard  sensors 
make  the  glider  heavier  and  require  onboard  decision  making  and  calculation  (Smith  et  al. 
2012),  the  glider  requires  more  energy  to  operate  and  consequently  its  mission  duration 
would  be  shortened.  Furthermore,  some  sensors  are  expensive;  for  example,  a  DVL  costs 
about  $50,000,  so  this  type  of  method  is  not  cost  effective.  Also,  since  DVL  estimates  the 
glider’s  position  by  integrating  its  estimated  velocity,  the  integrated  data  could  result  in 
cumulative  errors.  Singh  (2014)  attempted  to  improve  the  accuracy  of  the  Slocum’s 
underwater  position  by  developing  the  glider’s  dynamic  models  based  on  the  Extended 
Kalman  Filter  (EKF).  This  method,  however,  did  not  consider  the  ocean  currents  and 
glider’s  lateral  motion  as  well.  The  EKF  is  optimized  for  non-linear  process  and 
measurement  models.  In  addition,  its  prediction  operation  is  fast,  but  updating  the 
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measurement  is  slow  because  it  uses  matrix  inversion  (Pauli  et  al.  2014).  Thus,  the  error 
of  the  processing  delay  still  exists. 

The  various  filters,  including  the  EKF,  allow  the  gliders  to  estimate  its  position 
more  accurately  by  estimating  the  glider’s  state.  These  filters,  however,  also  have  some 
disadvantages.  The  KF  is  an  efficient  recursive  algorithm  that  can  predict  the  glider’s 
state  of  a  dynamic  system  with  measured  input  parameters  (Rajendra  and  Jannett  2007) 
such  as  pitch  angle  and  depth  by  the  onboard  sensors.  Since  the  KF  uses  the  measured 
input  parameters,  measurement  errors  by  onboard  sensors  can  occur,  and  the  errors  can 
be  accumulated  due  to  the  glider’s  long  mission  time.  Also,  the  KF  has  some  time  delays 
for  data  processing,  as  presented  by  Pauli  et  al.  (2014);  thus  the  delayed  data  causes 
errors  as  well.  Since  a  glider  can  only  update  parameters  consequent  to  its  position  by 
using  KF  when  it  locates  at  surface,  some  gliders  use  other  sensors  like  the  Inertial 
Navigation  System  (INS)  and  DVL  to  compensate  for  some  errors.  INS  uses 
accelerometers  and  gyroscopes  to  improve  DR  accuracy  by  integrating  measurements, 
and  DVL  can  measure  the  glider’s  velocity  relative  to  the  ground  (Woithe  et  al.  2011). 
Even  though  using  INS  and  DVL  is  more  accurate  than  only  using  the  KF  algorithm, 
calibration  and  initial  alignment  errors  can  still  exist  (Rajendra  and  Jannett  2007). 
Therefore,  it  is  essential  for  estimating  a  glider’s  precise  position  to  measure  input 
parameters  accurately  and  to  calibrate  onboard  sensors  precisely  when  additional  sensors 
are  used.  In  addition,  some  errors  can  be  reduced  by  making  the  glider  surface  frequently, 
but  this  would  negatively  impact  the  amount  of  ocean  data  collected  during  the 
underwater  mission. 

b.  Environmental  Model-Based  Methods 

The  environmental  model  focuses  on  the  effect  of  ocean  currents  for  estimating  a 
glider’s  position  by  combining  the  effects  of  currents  with  the  parameters  that  are  derived 
from  the  vehicle  model.  Since  the  glider  is  strongly  affected  by  ocean  currents,  the 
precision  of  its  position  can  be  increased  by  obtaining  more  accurate  underwater  currents 
infonnation. 
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Generally,  two  different  methods  can  be  used  to  obtain  the  ocean  current 
information.  One  is  measuring  the  in-situ  ocean  currents,  using  additional  onboard 
sensors  or  instruments  that  are  pre-installed  at  the  area  of  interest.  The  other  is  using  the 
ocean  model  data.  The  first  method  has  some  limitations  for  application  to  the  underwater 
gliders.  A  glider  can  be  equipped  with  additional  current  measurement  sensors  like  the 
Acoustic  Doppler  Currents  Profiler  (ADCP)  and  DVL  to  measure  the  in-situ  ocean 
currents  for  improving  the  accuracy  of  DR  algorithm.  DVL  can  calculate  the  ocean 
current  velocity  by  computing  the  velocity  of  the  ocean  water  and  the  seabed  relative  to 
the  glider.  ADCP  also  can  be  installed  within  the  Slocum  glider  to  estimate  real-time 
ocean  currents.  These  extra  onboard  sensors,  however,  have  some  drawbacks  as 
discussed  in  the  previous  section.  Moreover,  the  stability  of  ADCP  has  not  been  proven 
yet,  and  its  related  study  is  still  ongoing. 

Some  previous  studies  applied  the  ocean  model  prediction  data  to  plan  a  path  and 
to  design  trajectories  of  the  glider.  Smith  et  al.  (2010a,  b)  presented  the  effects  of  the 
ocean  currents  by  applying  the  ocean  model,  especially,  the  Regional  Ocean  Modeling 
System  (ROMS)  prediction  data  to  planning  a  path  and  designing  the  trajectory  of  the 
Slocum.  Their  simulations  show  that  there  was  approximately  a  50%  reduction  in  the 
errors  between  the  prescribed  waypoint  and  the  actual  surfacing  point  by  fusing  ROMS 
and  unscented  Kalman  Filter  (UKF).  As  presented  here,  the  previous  studies  considered 
the  ocean  currents  from  the  ocean  model  prediction,  but  these  studies  developed  the 
algorithm  in  planning  a  path  or  designing  trajectories  of  the  glider.  To  the  best  of  our 
knowledge,  no  previous  work  exists  for  localization  of  a  glider’s  underwater  position 
using  the  ocean  model  predictions.  Therefore,  this  thesis  attempts  to  develop  a  novel 
method  for  estimating  glider’s  underwater  position  using  the  ocean  currents  prediction 
data,  especially  from  the  Regional  NCOM  prediction. 

c.  Other  Methods 

Except  the  two  model-based  methods  discussed  in  the  previous  two  sections, 
some  researches  and  experiments  developed  other  methods  to  improve  the  accuracy  of  a 
glider’s  underwater  position.  But  these  studies  tried  to  estimate  the  position  by  using 
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additional  instruments  like  multiple  UUVs  and  acoustic  transponders  and  beacons  (Pauli 
et  al.  2014).  Somers  (2011)  studied  the  Doppler-based  localization  method  by  using 
multiple  Autonomous  Underwater  Vehicles  (AUVs),  and  Uffelen  et  al.  (2013)  studied  the 
method  using  broadband  acoustic  signals.  Even  though  these  methods  are  able  to  estimate 
a  glider’s  position  more  accurately,  they  are  not  cost-effective  methods  and  do  not 
consider  the  gliders’  dynamic  models.  Therefore,  the  details  of  these  methods  will  not  be 
covered  in  this  study. 

C.  CHAPTER  SUMMARY 

This  chapter  provided  motivations  for  developing  a  novel  method  for  estimating  a 
glider’s  optimal  underwater  trajectory  by  reviewing  general  characteristics  of  widely- 
used  gliders,  including  their  common  navigational  scheme  and  DR  algorithm,  each 
glider’s  navigational  features,  and  the  existing  model-based  methods. 

An  underwater  glider’s  position  accuracy  is  important  because  the  accuracy 
provides  spatiotemporally  higher  resolution  ocean  data  and  allows  the  glider  to  save 
energy  and  time,  as  well.  However,  since  the  glider  is  susceptible  to  the  ocean  currents 
and  lacks  inexpensive  and  efficient  localization  sensors,  estimating  its  precise  underwater 
position  is  difficult. 

The  three  widely-used  underwater  gliders,  the  Slocum,  the  Spray,  and  the 
Seaglider,  have  similar  appearances  and  a  common  navigational  scheme,  but  each  one 
has  unique  features  and  navigational  characteristics.  All  of  them  are  propelled  by  a 
buoyancy-driven  propulsion  system  and  are  designed  for  long  range  and  duration 
missions.  The  Spray,  the  Seaglider,  and  the  Slocum  Thennal  are  optimized  for  deep¬ 
water  missions.  Conversely,  the  Slocum  Electric  was  initially  optimized  for  shallow- 
water  coastal  missions,  but  recently  it  became  available  for  both  shallow  and  deep  water 
operations.  The  common  navigational  scheme  of  the  underwater  glider  is  mainly  reliant 
on  a  GPS  signal,  altimeter,  and  DR  algorithm.  The  DR  algorithm  is  the  most  frequently 
used  method  to  estimate  a  glider’s  underwater  position.  Even  though  the  algorithm 
provides  more  accurate  position  estimation  than  the  simple  linear  interpolation  method,  it 
has  some  margin  of  error.  The  main  error  is  caused  by  the  ocean  currents,  and  another 
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error  is  derived  from  the  DR  algorithm  itself  or  the  simplification  of  the  glider’s 
dynamics.  Also,  other  errors  come  from  the  parameters  measured  by  onboard  sensors  that 
relate  to  the  glider’s  dynamic  model.  The  Seaglider  uses  the  KF  algorithm  for 
compensating  its  navigation  error  by  estimating  the  velocity  of  the  ocean  currents.  The 
Slocum  and  the  Spray,  however,  do  not  estimate  the  ocean  currents,  so  their  estimated 
underwater  position  is  less  precise  than  that  of  the  Seaglider. 

Many  researches  and  studies  have  developed  model-based  methods  for  estimating 
a  glider’s  underwater  position  in  recent  decades.  These  methods  can  be  divided  into  two 
models:  the  vehicle  model  and  the  environmental  model.  Most  of  the  studies  focused  on 
the  vehicle  model-based  methods,  and  some  of  studies  developed  the  environmental 
model-based  methods  by  considering  the  ocean  currents  with  additional  onboard  sensors 
and  by  localizing  the  glider’s  position  with  additional  equipment.  These  developed 
methods  can  estimate  the  glider’s  underwater  position  more  accurately  than  the 
traditional  DR  algorithm,  but  still  have  some  limitations.  First  of  all,  most  studies  using 
the  glider’s  vehicle  model  have  limitations  in  explaining  the  effects  of  the  ocean  current 
because  the  model  does  not  consider  it.  Secondly,  the  methods  using  extra  onboard 
sensors  make  the  glider  heavier  and  require  more  energy  for  their  missions,  so  the 
glider’s  operation  duration  would  be  shortened.  Lastly,  the  methods  using  additional 
equipment  require  more  financial  backing,  and  the  equipment  is  sometimes  subjected  to 
the  complex  ocean  environment;  for  example,  the  ocean  currents,  bottom  topography. 
Therefore,  these  methods  are  less  cost-effective. 

The  Slocum  would  be  a  suitable  platform  for  littoral  operations  due  to  its  multi¬ 
depth  capability,  remarkable  versatility,  and  navigational  features  (Teledyne  Webb 
Research  2014).  Although  the  Slocum  has  some  advantages  for  littoral  operations,  it  still 
has  some  navigational  limitations.  Thus,  this  study  applies  the  Slocum  Electric  data  to 
develop  a  more  cost-effective  and  accurate  localization  method  based  on  the  ocean  model 
for  estimating  a  glider’s  optimal  underwater  trajectory.  This  novel  method  will  be 
discussed  in  Chapter  III. 
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III.  DATA  SET  AND  METHODOLOGY 


This  chapter  presents  the  methodology  to  develop  a  novel  approach  for  estimating 
a  glider’s  optimal  underwater  trajectory.  It  starts  with  a  description  of  the  Slocum  Electric 
data  set,  providing  general  information  on  field  experiments,  and  then  discusses  the 
Regional  NCOM  prediction  data  set.  Finally,  the  chapter  details  the  methodology. 

Subsurface  ocean  currents  must  be  considered  when  estimating  the  underwater 
position  of  a  glider.  Smith  et  al.  (2010a)  showed  the  effects  of  ocean  currents  in  planning 
a  path  of  the  Slocum  with  observation  of  approximately  a  50%  reduction  in  surfacing 
errors  by  incorporating  ROMS  predictions.  This  result  implies  that  the  ocean  model 
would  be  effective  in  estimating  a  glider’s  underwater  position.  If  real-time  ocean 
currents  information  is  available,  a  glider’s  underwater  position  could  be  estimated  more 
accurately,  but  it  is  difficult  to  obtain  the  infonnation  because  of  the  lack  of  in-situ 
measurement  sensors.  The  glider  must  be  equipped  with  extra  onboard  sensors  like 
ADCP  for  measuring  the  ocean  currents,  and  additional  instruments  must  be  pre-installed 
in  the  area  of  interest  to  measure  the  ocean  currents  or  to  estimate  the  glider’s  position  in 
near-real  time.  Before  installation  of  additional  instruments  for  ocean  current 
measurement,  we  propose  a  novel  method  for  estimating  a  glider’s  underwater  position 
by  incorporating  the  Regional  NCOM  prediction  data.  For  the  purpose  of  this  study,  it  is 
assumed  that  no  errors  exist  in  the  Regional  NCOM  data  and  the  measured  parameters 
from  onboard  sensors,  and  the  ocean  currents  structure  of  the  Regional  NCOM  is  similar 
with  that  of  in-situ  ocean  currents.  In  addition,  this  study  does  not  consider  the  dynamic 
or  mechanical  errors  of  the  glider. 

A.  DATA  SETS 

To  develop  the  novel  method,  two  data  sets  are  used:  the  Slocum  Electric  data  set 
and  the  Regional  NCOM  data  set.  The  Slocum  data  set  is  foundational,  providing  the 
trajectory  information  that  needs  to  be  corrected.  A  subset  of  the  Regional  NCOM  data, 
the  ocean  currents  prediction  data,  is  used  to  correct  the  inaccurate  trajectory. 
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1. 


The  Slocum  Electric  Data 


The  Slocum  Electric  (hereafter  referred  to  as  the  Slocum)  data  are  essential  for 
estimating  its  underwater  position.  The  data  are  obtained  from  a  Naval  Postgraduate 
School  class  field  experiment.  During  the  field  experiment  one  Slocum  was  deployed  at 
Monterey  Bay  (see  Figure  12)  for  three  days,  from  January  30,  2014,  to  February  1,  2014. 


Figure  12.  Experiment  Area  and  Trajectories  of  the  Slocum  (after  Google 

Earth). 


Actually,  this  field  experiment  was  not  intended  only  for  this  study,  but 
fortunately,  both  the  Slocum  and  Regional  NCOM  data  match  in  time  and  location,  so 
these  data  are  suitable  for  this  study.  During  the  three  days  of  deployment  time,  the 
Slocum  performed  68  dives.  The  average  straight  distance  of  DR  trajectory  was  about 
362  meters,  with  a  minimum  of  112  meters  and  maximum  of  623  meters,  depending  on 
the  prescribed  depth  and  location.  The  average  time  of  a  single  dive  was  about  20.5 
minutes,  and  the  average  time  of  staying  at  surface  was  about  10  to  15  minutes.  In 
addition,  the  average  of  maximum  depth  was  about  95  meters,  with  maximum  depth  of 
about  198  meters.  After  finishing  the  three-day  mission,  the  Slocum’s  basic  navigation 
data  were  processed  with  MATFAB.  The  data  included  time;  pressure;  CTD;  and  the  data 
for  three  main  positions,  i.e.,  GPS  fixed  position,  DR  position,  and  linear  interpolation 
position.  To  estimate  the  glider’s  optimal  underwater  trajectory,  we  assume  that  there  are 
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no  errors  in  DR  trajectory.  Also,  we  do  not  consider  surface  drift,  so  we  only  attempt  to 
estimate  the  underwater  trajectory  for  each  dive. 

2.  The  Regional  NCOM  Data 

NCOM  is  a  free-surface  general  circulation  model  (Ngodock  and  Carrier  2014) 
based  primarily  on  two  existing  ocean  circulation  models,  the  Princeton  Ocean  Model 
(POM)  (Martin  et  al.  2009)  and  the  Sigma/Z-level  Model  (SZM)  (Martin  et  ah,  1998).  It 
is  also  based  on  primitive-equation  and  using  the  hydrostatic,  Boussinesq,  and 
incompressible  approximations  (Martin  et  al.  2009).  NCOM  is  a  real-time  data- 
assimilating  global  ocean  nowcast  and  forecast  system.  It  has  been  developed  at  the 
Naval  Research  Laboratory  (NRL)  and  transitioned  to  the  NAVOCEANO  to  support 
naval  operations,  operational  activities,  and  other  researches  (Rhodes  et  al.  2002).  It  has 
been  used  in  global-  and  basin-scale  circulation  applications  (Barron  et  al.  2003,  2004; 
Kara  et  al.  2006). 

NCOM  has  two  types  of  model:  the  Global  and  the  Regional  model.  The  Global 
NCOM  (GNCOM),  however,  has  been  replaced  by  the  operational  Global  HYbrid 
Coordinate  Ocean  Model  (HYCOM)  (Ocean  Prediction  Center  (OPC)  2014).  HYCOM 
has  a  horizontal  resolution  of  1/12  degree  (about  9  kilometers)  (Metzger  et  al.  2014)  and 
40  depth  levels,  and  it  uses  hybrid  (isopycnal/sigma/z-level)  coordinates  in  the  vertical 
(National  Oceanic  and  Atmospheric  Administration  (NOAA)  2014).  The  HYCOM  data 
consist  of  sea  surface  temperature  (SST),  sea  surface  height  (SSH),  eastward  and 
northward  currents  and  subsurface  temperature  and  salinity  (NOAA  2014;  NRL  2014). 
HYCOM  assimilates  data  from  satellite  and  in-situ  observation,  using  the  Navy  Coupled 
Ocean  Data  Assimilation  (NCODA)  system  (NRL  2014).  HYCOM  provides  boundary 
conditions  to  the  Regional  NCOM  (NOAA  2014).  The  Regional  NCOM  includes  the 
U.S.  East  coast,  the  Southern  California  coast,  Hawaiian  coasts,  and  the  Gulf  of  Mexico 
and  Caribbean  Seas  (OPC  2014). 

In  this  thesis,  the  Regional  NCOM  data  of  Southern  California  coast  are  used. 
This  data  set  has  a  high  resolution  of  1/30  degree,  about  3.6  kilometers.  Every  day,  the 
U.S.  Navy  produces  a  seven-day  forecast  (Metzger  et  al.  2014)  of  the  Regional  NCOM. 
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The  data  is  updated  at  OOZ  with  the  four-day  (96  hours)  forecast  and  three-hour 
increments  (NOAA  2014).  The  Regional  NCOM  has  the  grid  data  which  consist  of  many 
variables,  especially  the  time;  position  (longitude  and  latitude);  depth;  and  the  eastward 
( u ncom  )  ar|d  northward  ( vncom  )  component  of  ocean  current  velocity.  Even  though  the 
Regional  NCOM  has  high  resolution  data,  it  requires  being  reprocessed  to  match  with  the 
Slocum  data.  After  the  processing  is  done,  these  data  can  be  applied  to  our  novel  method. 
Details  of  data  processing  of  the  Regional  NCOM  are  described  in  the  following  section. 

a.  Data  Processing  of  the  Regional  NCOM 

The  Slocum  data  is  point  data,  while  the  Regional  NCOM  data  is  grid  data  (see 
Figure  13),  so  data  processing  is  required  to  match  both  data  sets  with  time  and  location. 


Depth: 
40  levels 


Lonaitude:  424  points 


Figure  13.  Schema  of  Grid  Data  of  the  Regional  NCOM  for  One  Step  Time 

Interval. 


The  water  velocity  of  Regional  NCOM  is  divided  into  eastward  ( uncom  )  and 
northward  ( vncom  )  components,  and  these  components  depend  on  the  4-D  variables: 

longitude,  latitude,  depth,  and  time.  The  time  duration  of  prediction  is  three  hours  in  one 
data  file.  Given  every  three-hour  time  interval,  the  longitude  is  about  14  degrees  long 
(from  1 10°54’  W  to  125°00’  W)  dividing  linearly  by  424  points,  and  the  latitude  is  about 
15  degrees  long  (from  25°00’  N  to  40°06’  N)  dividing  linearly  by  454  points,  which 
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means  the  resolution  of  position  data  is  1/30  degree,  about  3.6  kilometers.  The  depth  data 
is  divided  non-linearly  by  40  levels  from  surface  to  5,000  meters.  The  first  20  steps  are 
from  surface  to  100  meters,  and  then  the  depth  is  rapidly  increasing  to  maximum  depth. 
Given  this  original  Regional  NCOM  data,  we  use  linear  interpolation  for  the  data 
processing  to  match  both  the  Slocum  and  the  Regional  NCOM  data, 


m= 
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Now,  an  example  of  the  method  for  calculating  u  at  specific  time  is  presented. 
Let  us  compute  the  data  of  u  at  time  10:00  with  the  same  position  and  depth.  Since  the 
time  interval  is  three  hours,  choose  two  boundary  times,  9:00  and  12:00,  which  include 
10:00.  If  the  velocity  at  9:00  and  12:00  are  0.5  ms'1  and  0.3  ms'1,  respectively,  the 
velocity  at  10:00  can  be  calculated  as  follows. 


«(10)  - 


0.5 


12-10 

12-9 


+ 


0.3 


10-9 

12-9 


=  0.5— +  0.3  —  «  0.4333  ms~l 
3  3 


This  linear  interpolation  method  is  applied  to  all  data  sets. 


B.  METHODOLOGY 

In  this  section,  a  novel  method  will  be  developed  through  two  phases.  Phase  One 
proves  the  importance  of  the  ocean  currents  by  comparing  distance  c/,  to  d2  .  <7,  is  the 

distance  between  the  actual  resurfaced  position  (Pgps)  and  final  estimated  position  (Pd,)  of 
DR  trajectories  (T i),  and  d2  is  the  distance  between  Pgps  and  the  final  estimated  position 
(P ncom)  of  updated  DR  trajectories  (T2).  We  also  refer  to  this  trajectory  as  DRNCom- 
DRncom  trajectory  is  obtained  from  Ti  by  combining  its  DR  position  data  and  the 
Regional  NCOM  data  with  the  simple  composite  trapezoidal  rule,  which  is  a  type  of 
numerical  integration. 

Phase  Two  develops  a  novel  correction  method  to  estimate  the  glider’s  optimal 
underwater  trajectory  by  correcting  with  the  depth-averaged  and  depth-dependent 
correction  methods  and  by  using  the  simple  composite  trapezoidal  rule.  Both  correction 
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methods  would  increase  precision  for  estimating  the  glider’s  trajectory,  but  to  develop  a 
better  method,  they  will  be  compared.  The  depth-averaged  correction  method  idea  is 
originated  from  the  depth-average  currents  calculation  method.  The  depth-averaged 
correction  method  calculates  the  water  displacement  by  dividing  the  total  distance  of 
error  (AXW  ,AYW)  by  time  ratio  ( tn/tN  )  between  total  transit  time  ( tN  )  and  at  each  time 
(  tn  ).  Since  this  method  corrects  the  trajectory  by  linearly  calculating  the  water 

displacement,  we  refer  to  it  as  the  Linear  Correction  Method  (LCM).  The  depth- 
dependent  correction  method  is  a  novel  approach  for  estimating  a  glider’s  trajectory, 
correcting  the  trajectory  by  using  the  spatiotemporally  dependent  ocean  prediction  data. 
Unlike  the  depth-averaged  correction  method,  this  method  calculates  the  water 
displacement  by  using  the  ocean  currents  prediction  data  that  depend  on  the  depth  and 
time  consequent  to  the  location,  so  the  solution  cannot  be  calculated  linearly.  Thus,  we 
referred  to  it  as  the  Non-linear  Correction  Method  (NCM). 

For  this  study,  the  Phase  Two  is  divided  into  four  steps.  In  Step  One  and  Two,  the 
corrected  trajectories  by  two  correction  methods  are  described.  In  Step  Three,  an  optimal 
trajectory  is  estimated,  and  the  corrected  trajectories  are  compared  to  the  optimal 
trajectory  in  Step  Four.  Since  this  study  focuses  only  on  the  impact  of  the  ocean  currents 
instead  of  other  factors  like  the  glider’s  dynamic,  mechanics,  and  measurements,  these 
errors  are  neglected. 

1.  Phase  One:  Proof  of  the  Ocean  Effects  on  Localization 

The  purpose  of  this  section  is  to  prove  the  hypothesis  that  the  underwater  position 
of  a  glider  is  strongly  affected  by  the  ocean  currents,  and  its  position  estimation  can  be 
improved  by  applying  the  ocean  model  prediction  data.  Once  this  hypothesis  is  proven, 
the  Regional  NCOM  prediction  data  will  be  applied  to  estimating  the  glider’s  optimal 
underwater  trajectory  in  Section  2. 

To  verify  the  hypothesis,  we  compare  the  distance  between  the  resurfaced  and 
estimated  final  positions  {Pdr,  Pncom,  and  Pgps)  of  each  trajectory  (Ti,  T2,  and  T3)  in  each 
dive.  These  trajectories  are  defined  as  follows:  The  first  trajectory  (T i)  is  calculated  from 

the  traditional  DR  algorithm  of  the  glider  as  we  described  in  Chapter  II.  The  second 
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trajectory  (T2)  is  updated  by  combining  the  glider’s  position  data  on  Ti  with  the  water 
displacement  that  is  calculated  from  the  Regional  NCOM  data  on  T3  by  numerical 
integration.  The  third  trajectory  (T3)  is  calculated  by  linear  interpolation  between  two 
GPS  fixed  positions  in  each  dive.  Though  the  linear  interpolation  is  not  realistic,  this 
trajectory  is  used  to  estimate  T2.  Since  the  glider’s  actual  underwater  position  cannot  be 
known,  we  assume  the  glider  follows  T3.  In  addition,  since  the  time  of  each  point  on  Ti 
and  T3  can  be  matched,  the  water  velocity  of  the  Regional  NCOM  data  at  each  position 
on  T3  also  can  be  used  to  estimate  the  water  displacement  consequent  the  position  of  T2. 
Ti  and  T3  are  calculated  from  the  glider’s  algorithm,  so  they  are  set  as  reference 
trajectories  to  estimate  the  T2.  To  estimate  the  position  of  T2,  both  the  glider’s 
displacement  and  the  water  displacement  are  required  to  be  calculated.  The  glider’s 
displacement  was  already  calculated  from  the  DR  algorithm,  but  the  water  displacement 
is  required  to  be  calculated.  Since  the  Regional  NCOM  prediction  data  are  grid  data,  the 
linear  interpolation  is  used  to  match  the  time  and  depth  with  each  point  on  T]  as 
described  in  Section  A.  Pgps  on  T3  is  merely  obtained  from  the  GPS  fixed  information. 
With  these  three  resurfaced  positions,  we  compared  the  distance  between  Pdr  and  Pgps  ( d1  ) 

to  Pncom  and  Pgps  ( P  )  • 

Generally,  the  displacement  ( Rn )  of  marine  vehicles  can  be  simply  calculated  by 
integration  with  their  velocity  ( U(t) )  during  the  transit  time  ( tn ). 

R„  =  f  U(t)  dt,  n  =  \,2,--,N  -\,N  (11) 

Jt0 

Since  the  glider’s  underwater  velocity  (Ug(t))  is  relative  to  the  ocean  currents’ 
velocity  ( U ncom(t) ),  their  total  underwater  velocity  (  Utot (/) )  can  be  expressed  as 

u(0  =  uto((0  =  ug(0 + u„com(0  (12) 

By  substituting  U,o,(0  in  Equation  (12)  into  Equation  (11),  the  glider’s 
underwater  displacement  can  be  estimated  analytically  by  numerical  integration 

R,  =  !>„«)*  =  f''ui(r)*  +  pura(0*  (13) 

J  *0  Jt0  Jt0 
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In  this  study,  however,  displacement  of  the  glider  can  be  estimated  without 
calculation  because  its  displacement  ( f  U(f)  dt)  is  already  calculated  as  longitude  and 

Jt0 

latitude,  and  thus,  the  only  thing  required  to  calculate  is  the  water  displacement 
(  f  U „com(t)dt).  In  this  study,  the  analytical  solution  can  be  simply  calculated  as  a 
numerical  solution  by  using  the  composite  trapezoidal  rule, 

fvU  (t)dt=['  U  (t)dt  +  f2U  (t)dt  +  —  +  f'w  U  (t)dt 

J  ncom  \  /  J  ncom+J  J  ncom  \  /  J  ncom\  s 

•I  t0  •>  [q  J  l j  •/  tN_\ 

=  |(U0+U,)  +  ^(Ul+U,)  +  ...  +  |(U„_1+U„)  (14) 

-~(U0  +2U,  +2U,  H  +2U^_,  +UWJ,  h  =  ‘-E-I± 

This  rule  can  be  used  when  the  value  of  each  time  interval  (At)  is  constant,  but 
At  in  this  study  is  not  constant.  However,  the  mean  values  of  the  time  interval  and  the 
velocity  difference  between  each  point  are  small:  AT  =  4.346sec,  u  =2.06x\0  6 ms^1 , 
v  =-1.55x10  W1,  respectively,  so  the  error  would  be  small.  Therefore,  the  simple 
composite  trapezoidal  rule  can  be  applied  to  calculate  the  water  displacement  expressed 
as 


f '  U „com ( t )  dt  =  f '  Uncom (t)dt +  ■■■+[ '  Uncom (t)  dt 

"  'o  *  ‘0  "  lN-\ 

=  ^L(U„  +  U1)  +  -  +  ^(U„_,+U„),  At,  2, -,N) 


(15) 


Generally,  a  glider’s  two-dimensional  (2-D)  movement  consists  of  eastward  (x) 
and  northward  ( y  )  components,  so  does  the  water  displacement.  Also,  water  velocity 
( U  )  consists  of  eastward  ( u  )  and  northward  ( v  )  components  corresponding  to 

x  ncom  '  x  ncom  '  x  r  ncom  /  r  i  o 


the  water  displacement  ( Sx.  ,  8y'-  )  at  each  time  interval  ( At. ), 

cli  At 

S*7  =  I ,  Uncom  (0  dt  =  — ^  (U[  ,  +11,  ),  At ,  =  t,  ~  tl  A 

J  +i  2 

rt.  At 

Sy'i  =  j,  V„com  (0  dt  —  — L  (vM  +  vf ),  A t:  =  t,  -  G, 


(16) 
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Then,  the  water  displacement  at  time  tn  will  be 


a *;=2>; 

i= 1 


(17) 


Since  the  unit  of  calculated  water  displacement  in  Equation  (17)  is  meters,  it  is 
required  to  be  converted  into  units  of  degree  by  using  Equation  (18)  to  estimate  the 
underwater  position  of  To.  Note  that  we  assume  the  one  degree  in  both  longitude  and 
latitude  as  110  kilometers,  i.e.,  l.lxlO5  meters. 

/on"  =  Ax"  /  (1 . 1  x  1 05 ) 

/a/"  =  Ay"  /  (l.lxlO5) 


(18) 


Then,  each  position  of  To  can  be  obtained  by  adding  the  water  displacement  to  the 
longitude  and  latitude  of  the  glider  displacement.  Finally,  the  glider’s  underwater  position, 
which  is  updated  by  the  Regional  NCOM  prediction,  can  be  obtained  as 


lonnncom  =  londnr  +  Ion ",  n  =  1,  2,  •  •  -,  N  - 1,  N 

latnncom  =  latdnr  +  latZ ,  n  =  1,  2,  •  •  ■ ,  N  - 1,  N 


(19) 


where  Ion'll ,  lonnncom  and  lath  ■>  are  the  longitude  and  latitude  of  T i  and  T2  at  each 


time  ( tn  )■  To  calculate  the  distance  dx  and  d2 ,  first  calculate  the  distance  with  longitude 
and  latitude  by  using  the  distance  formula;  then  convert  the  unit  of  distance  from  degrees 
to  meters. 


Let  the  longitude  and  latitude  of  Pgps,  Pdr,  and  Pncom  be  longps ;  londr ;  lonncom  ; 
lat?ps ;  latdr ;  and  latncom  ,  respectively.  Then  each  distance  dx  and  d2  is  calculated  as 

d\  =  J(/ongPs  ~l°ndrf  +  (latgps  ~latdrf  (1-lxlO5) 

V _ _ _  _ _  (20) 

d2  ^  -l0nncom)2  +  (lat gPs  ~latncomf  (hlxlO') 

After  calculating  the  distance  d}  and  d 2 ,  compare  the  result  to  prove  how  well 
estimated  the  trajectory  is  by  considering  the  ocean  model  prediction  data. 
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2.  Phase  Two:  Trajectory  Correction  Method  Development 

To  develop  a  novel  method  for  estimating  a  glider’s  optimal  underwater  trajectory, 
two  correction  methods  are  developed:  the  Depth-averaged  Correction  Method,  which  is 
referred  to  as  the  Linear  Correction  Method  (LCM)  and  the  Depth-dependent  Correction 
Method,  which  is  referred  to  as  the  Non-linear  Correction  Method  (NCM).  Newly 
corrected  trajectories  would  be  produced  by  combining  DR  trajectory  with  the  ocean 
model  prediction  data.  These  corrected  trajectories  will  be  developed  in  the  first  two  steps, 
and  an  optimal  trajectory  will  be  estimated  in  Step  Three.  Then  two  corrected  trajectories 
will  be  compared  to  the  optimal  trajectory  in  Step  Four  to  detennine  how  well  the 
methods  are  developed. 

a.  Depth-Averaged  Correction  Method:  Linear  Correction  Method  (LCM) 

The  depth-averaged  correction  is  not  a  newly  developed  method,  but  the 
application  of  this  method  to  estimate  a  glider’s  underwater  trajectory  would  be  a  new 
attempt.  This  method’s  inspiration  is  derived  from  the  depth-averaged  currents 
calculation  method.  Thus,  this  method  is  similar  to  the  calculation  method  of  the  depth- 
averaged  currents.  Figure  3  promotes  a  better  understanding  of  this  method. 

When  the  glider  resurfaces  after  each  dive,  Pgps  is  the  only  exact  position,  and  Pd,- 
is  the  estimated  position  by  the  DR  algorithm.  Thus,  a  position  error  between  Pgps  and  Re¬ 
occurs.  Since  the  DR  algorithm  uses  only  the  glider’s  velocity,  which  is  relative  to  the 

water  velocity,  and  does  not  consider  the  ocean  currents,  the  total  distance  ( AR")  of  the 
error  accounts  for  the  ocean  current  effects,  i.e.,  the  water  displacement.  Therefore,  more 
precise  estimation  of  the  glider’s  underwater  trajectory  would  be  obtained  by  adding  the 
ocean  currents  displacement  to  the  DR  position.  The  total  water  displacement  consists  of 

eastward  (AX'1  )  and  northward  (AyM)  components,  and  the  glider’s  total  transition  time 
is  ( tN ).  Before  detailing  the  method,  let  us  assume  that  the  glider’s  actual  underwater 
trajectory  follows  on  the  linear  interpolation  trajectory  (T3)  because  no  one  knows  the 
actual  underwater  trajectory  of  the  glider.  With  this  assumption,  it  is  possible  to  calculate 
the  water  displacement  ( AxJ  ,  A y”  )  at  each  time  ( tn  )  by  using  the  depth-averaged 
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correction  method.  The  water  displacement  at  each  time  can  be  simply  calculated  from 
the  total  water  displacement  (AX’1  ,AyH)  divided  by  the  time  ratio  between  the  total 
transit  time  and  at  each  time  ( tn  /tN  ).  Then,  the  water  displacement  at  each  time  can  be 

expressed  linearly  as 

Ax”’  =  AXW  —  , 

H  t 

lN 

Ay”  =  AYW  —  , 
t N 

Since  the  water  displacement  is  calculated  linearly  consequent  to  the  estimated 
position,  we  refer  to  this  method  as  the  Linear  Correction  Method.  As  in  the  first  phase, 
since  the  unit  of  each  DR  position  data  ( Ion f  ,  latf  )  is  degrees  and  the  unit  of  water 


w  =  l,  2, ---,  N-\,N 
n  =  \,  2,--,N-\,N 


(21) 


displacement  is  meters,  the  latter  needs  to  be  converted  to  degrees  ( Ion™' ,  latcnor )  by 

/o<°'=Ax;/(l.lxl05) 
latc°r  =  Ay;  / (l.lxlO5) 

Then,  the  estimated  underwater  position  of  the  glider  ( lonenst ,  lat^' )  at  each  time 
can  be  expressed  as 

lonesl  =londr  +loncor 

"  "  n  (23) 

7  ^est  7  idr  .  7  j.cor 

tatn  =  latn  +  Latn 

We  can  estimate  the  glider’s  underwater  position  using  this  simple  and  fast 
correction  method.  This  method,  however,  does  not  consider  the  actual  ocean  currents 
since  the  real  position  cannot  be  obtained.  Moreover,  it  calculates  the  water  displacement 
linearly  and  does  not  consider  the  ocean  currents’  dependence  on  depth  and  time. 
Therefore,  the  more  precise  method  using  depth-dependent  correction  will  be  discussed 
in  the  following  section. 


b.  Depth-Dependent  Correction  Method:  Non-linear  Correction  Method 
(NCM) 

Depth-dependent  correction  is  a  novel  method  which  considers  the  spatiotemporal 
effects  of  the  ocean  current,  so  it  is  expected  that  a  more  precise  underwater  position  of 
the  glider  would  be  estimated  with  this  method.  Like  the  depth-averaged  correction 
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method,  the  depth-dependent  correction  method  considers  the  effects  of  ocean  current, 
but  the  latter  uses  the  spatiotemporally  dependent  data  of  ocean  model  prediction  to 
estimate  more  precise  water  displacement  rather  than  the  depth-averaged  water 
displacement  for  estimating  a  glider’s  underwater  trajectory.  Also,  the  assumption  that 
the  glider’s  actual  underwater  trajectory  follows  on  the  linear  interpolation  trajectory  is 
still  valid.  The  primary  difference  between  the  two  methods  is  the  calculation  method  of 
the  water  displacement.  The  depth-averaged  correction  method  calculates  the  water 

displacement  linearly  by  multiplying  total  displacement  (AAM',A7tv)  by  time  ratio  (tn/tN  ) 
at  each  time  ( tn )  as  shown  in  Equation  (21).  The  depth-dependent  correction  method,  on 
the  other  hand,  calculates  the  water  displacement  by  using  the  Regional  NCOM 
prediction  data,  which  is  dependent  on  the  depth  of  the  glider  (  z  )  at  each  time- 
consequent  location.  Since  this  method  calculates  the  water  displacement 
spatiotemporally,  the  solution  is  not  linear.  Therefore,  we  refer  to  this  method  as  the  Non¬ 
linear  Correction  Method  (NCM). 

In  this  method,  Equation  (11)  to  Equation  (13)  can  be  applied  to  estimate  the 
glider’s  underwater  position,  and  the  glider’s  subsurface  displacement  is  also  calculated 
by  the  DR  algorithm  with  the  Slocum  data.  Thus,  the  only  aspect  that  needs  to  be 
calculated  is  the  water  displacement.  The  water  displacement  ( Ar“  )  at  time  ( tn  )  is 
calculated  from  the  Regional  NCOM  data  on  T3  by 

AC  =  f  U„com (0  dt,  n  =  l,2,-,N-l,N  (24) 

The  Regional  NCOM  prediction  data  have  4-D  variables,  i.e.,  U ncom(x,  y,  z,t) , 
and  the  water  velocity  consists  of  the  eastward  (  u  )  and  the  northward  ( v  ) 

component,  as  already  mentioned  in  Section  1  in  this  chapter.  In  this  section,  however, 
we  will  only  present  the  method  to  calculate  eastward  displacement.  The  northward 
displacement  can  be  calculated  by  substituting  v  with  u  in  each  equation.  The 

water  displacement  can  be  simply  calculated  by  the  simple  composite  trapezoidal  rule  of 
numerical  integration,  which  is  the  same  as  in  Equation  (16)  and  Equation  (17). 
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Therefore,  the  water  displacement  of  the  eastward  component  ( Ax’’  )  at  time  ( tn  )  is 
calculated  as 


Ax”’  =  u  (t)  dt-  T 

n  I  ncom  V  '  /  j 

'°  i= i 


At,.  . 

—  (wM+n,) 


n 

Z(  ncom  \  j  \ 
i= 1 


where  w 


ncom  _  R  1  "t"  R 


and  At,.  =  t,.  - 1,  , 


(25) 


In  Equation  (25),  however,  the  depth-dependent  water  velocity  of  the  Regional 
NCOM  prediction  data  is  required  to  calculate  the  water  displacement,  using  the  depth- 
dependent  correction  method.  In  addition,  the  water  velocity  of  the  Regional  NCOM 
prediction  data  depends  on  the  time-consequent  location,  so  the  eastward  component  of 
the  water  velocity  can  be  expressed  as  u"com{zi,ti) .  Even  though  the  depth-dependent 


water  velocity  is  used  to  correct  the  trajectory,  the  corrected  trajectory  would  not  follow 
the  real  trajectory,  and  the  final  position  of  the  corrected  trajectory  would  not  match  with 
Pgps,  like  To  in  Figure  2.  Therefore,  some  equations  for  correcting  the  water  velocity  will 
be  discussed  to  match  the  final  position,  and  thus  to  correct  the  trajectory  more  precisely. 


The  eastward  component  of  the  water  velocity  of  the  Regional  NCOM  prediction 
data  at  each  time  can  be  expressed  in  three  parts  as 


„  .  ncom  /  ,  \ 

ui  (zp  6)  = 


ncom  /  ,  \  ncom  /a  ^  \ 

«,  (zp  0 u,  (°c,). 


(26) 


ncom  /  a  ^  \  llidA 

»,  (0,*,)  wmax 

where  umax  =  max  [u"com  (0,  t{ )] ,  i  =  1,  2,  •  •  • ,  N 

The  first  part  of  the  right-hand  side  shows  the  velocity  ratio  between  the  surface 
and  the  glider’s  depth  at  time  tt ,  and  the  second  part  presents  the  surface  velocity  ratio 

between  the  position  of  maximum  surface  velocity  during  transit  time  of  one  dive  ( tN  ) 


and  the  position  of  the  glider  at  time  .  We  refer  to  these  two  parts  as  z  factor  ( f: ),  and  t 
factor  ( f-  ),  respectively. 


fr  _  uTM  t  =  ur; 

J  i  ncom  /r\  ,  \  ?  J  n 


'(0,6) 


U; 


'(0,6) 


then,  u"com(zi,  t.)  =  fz  f„‘un 


(27) 
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The  third  part  indicates  the  maximum  surface  velocity  during  the  transit  time  of 
one  dive,  and  the  maximum  velocity  has  a  constant  value.  If  the  value  of  maximum 
velocity  is  the  same  as  shown  Equation  (26),  the  corrected  trajectory  would  be  exactly 
same  with  the  results  of  the  Phase  One.  Thus,  the  velocity  of  the  Regional  NCOM  data 
for  trajectory  correction  will  be  redefined  more  realistically  as 


urM=ftfK 


cor 

max 


(28) 


Moreover,  since  no  one  knows  the  actual  trajectory  of  the  glider  consequent  the 
value  of  the  maximum  surface  velocity,  a  new  equation  that  does  not  use  the  maximum 
surface  water  velocity  is  required.  For  this  study,  it  is  assumed  that  there  is  no  error 
between  the  ocean  current  structure  of  in-situ  data  and  that  of  the  Regional  NCOM 
prediction  data.  Given  this  assumption,  we  can  apply  z  factor  and  t  factor  to  correct  the 
trajectory.  Let  the  total  eastward  water  displacement  at  final  position  ( tN  )  during  one 


dive  be  AX" .  Substitution  of  Equation  (28)  into  Equation  (25)  leads  to 

N  (  N  "\ 


Ax-=y„r(z„»,)Af,= 


u. 


(29) 


i=l  V  i= 1  J 

After  that,  the  maximum  surface  velocity  to  correct  the  trajectory  is  re-expressed 


as 


.  cor 

u  „  =  - 


AT 


max  N 


(30) 


Zf'fi  A‘i 


i= 1 


Similarly,  the  water  displacement  of  each  trajectory  at  time  tn  is  given  by 


V  (-1  J 


i= 1 


(31) 


Finally,  the  water  displacement  of  each  trajectory  at  time  tn  can  be  re-expressed 


as 


A*:  =  AX" 


A», 

/•=! _ 

Af, 


1=1 


(32) 
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After  calculating  the  corrected  water  displacement,  processing  for  estimating  the 
glider’s  underwater  trajectory  is  the  same  with  the  LCM  as  shown  in  Equation  (22)  and 
Equation  (23). 


c.  Application  to  Four  Steps 

This  section  describes  the  application  of  two  developed  correction  methods  to  DR 
trajectory  correction.  In  first  two  steps,  DR  trajectory  is  corrected  by  using  two  correction 
methods  without  comparing  because  there  is  no  reference  trajectory  to  determine 
improvement  of  estimated  trajectory.  In  Step  Three,  an  optimal  trajectory  is  estimated  as 
a  reference  trajectory  by  using  iteration  method,  incorporating  the  depth-dependent 
correction  method.  Then,  two  corrected  trajectories  are  compared  to  a  reference  trajectory 
in  Step  Four.  By  comparing  these  three  trajectories,  it  would  be  possible  to  determine 
which  correction  method  is  better  for  estimating  a  glider’s  optimal  underwater  trajectory. 
The  brief  explanation  of  each  step  is  as  follows. 

(1)  Step  One:  The  Depth-averaged  Correction  with  DR 

In  this  step,  DR  trajectory  (T i)  is  corrected  using  the  depth-averaged  correction 
method.  The  schema  of  this  method  is  the  same  as  shown  in  Figure  3.  This  corrected 
trajectory  is  referred  to  as  Ci.  The  Ci  is  expected  to  be  a  less  accurate  trajectory  than 
other  corrected  trajectory. 

(2)  Step  Two:  Depth-dependent  Correction  with  DR 

In  this  step,  Ti  is  corrected  using  the  depth-dependent  correction  method.  The 
schema  of  this  method  is  the  same  as  shown  in  Figure  4.  This  corrected  trajectory  is 
referred  to  as  C2.  C2  is  expected  to  be  a  more  accurate  trajectory  than  Ci. 

(3)  Step  Three:  Optimal  Trajectory  Development  using  Iteration  Method  with 
the  C2 

In  this  step,  more  precise  trajectory  of  the  glider  is  estimated  by  iterating  the  C2 
until  relative  error  converges  to  10-5.  This  corrected  trajectory  is  referred  to  as  C3.  C3  is 
expected  to  be  the  most  precise  trajectory  among  the  three  trajectories.  Thus,  it  is  also 
referred  to  as  an  optimal  trajectory. 
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(4)  Step  Four:  Comparison  Three  corrected  Trajectories 

In  this  step,  two  corrected  trajectories  will  be  compared  to  an  optimal  trajectory 
(Ci  and  C3,  C2  and  C3)  to  determine  which  correction  method  is  more  precise  for 
estimating  a  glider’s  optimal  underwater  trajectory. 

C.  CHAPTER  SUMMARY 

To  increase  the  precision  of  a  glider’s  trajectory,  a  novel  method  was  developed 
by  applying  the  Regional  NCOM  data  through  two  phases  and  four  subset  steps  in  Phase 
Two. 

In  Phase  One,  we  attempted  to  show  the  influence  of  ocean  currents  by  applying 
the  Regional  NCOM  data  to  the  DR  trajectory.  In  Phase  Two,  we  aimed  to  develop  two 
correction  methods:  the  depth-averaged  and  the  depth-dependent  correction  methods  for 
estimating  the  glider’s  optimal  underwater  trajectory.  Even  though  both  correction 
methods  consider  the  ocean  currents,  the  former  corrects  the  glider’s  trajectory  linearly; 
and  this  linear  correction  method  would  be  relatively  less  accurate.  Thus,  to  correct 
trajectory  more  precisely,  we  developed  the  latter  method.  Since  the  latter  method 
considers  the  ocean  currents  dependence  on  depth  and  time-consequent  location,  it  could 
estimate  a  glider’s  underwater  trajectory  more  accurately  than  the  former  one.  With  these 
two  developed  methods,  we  corrected  the  DR  trajectory  by  applying  them  to  two 
different  steps.  Then,  two  correction  methods  are  evaluated  by  comparing  two  corrected 
trajectories  to  an  optimal  trajectory.  In  Chapter  IV,  the  results  of  these  correction 
methods  will  be  analyzed  by  comparing  the  corrected  trajectories. 
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IV.  RESULTS  AND  DISCUSSION 


This  chapter  provides  results  of  this  study,  analyzing  corrected  trajectories  by 
applying  the  two  correction  methods  that  were  described  in  Chapter  III.  Detailed  results 
of  each  phase  and  step  will  be  discussed,  as  well  as  the  relevance  of  this  study  to  the 
ROKN. 


A.  RESULTS 

The  analysis  is  divided  into  two  subsections.  In  the  first  section,  the  DR  (T i)  and 
the  DRNcom  (T2)  trajectories  are  compared  to  each  other  to  show  that  the  Regional 
NCOM  prediction  data  have  an  impact  on  improving  the  accuracy  of  a  glider’s  DR 
trajectory.  After  obtaining  updated  trajectories,  the  distances  between  final  points  of  each 
trajectory  Pgps  and  P  dr  ( dx  ),  Pgps  and  PnCom  ( d2 )  are  compared  to  each  other.  In  the  second 
section,  it  is  determined  which  method  is  a  better  correction  method  for  estimating  a 
glider’s  optimal  underwater  trajectory  throughout  the  four  steps.  Step  One  deals  with 
results  of  DR  trajectories,  which  are  corrected  by  the  depth-averaged  correction  method 
(Ci).  Step  Two  provides  results  of  DR  trajectories,  which  are  corrected  by  the  depth- 
dependent  correction  method  (C2).  Since  there  is  no  actual  trajectory  that  can  compare 
with  the  Ci  and  C2  to  prove  their  estimation  precision,  an  optimized  trajectory  (C3)  is 
estimated.  After  obtaining  the  C3,  the  Ci  and  the  C2  are  compared  to  C3  to  determine 
which  method  is  a  better. 

1.  Result  of  Phase  One 

This  section  presents  the  importance  of  the  ocean  current  effects  in  estimating  a 
glider’s  underwater  trajectory  by  comparing  the  distance  between  dx  and  d2  for  each 
dive.  Figure  14  shows  an  example  of  the  results,  and  more  figures  of  results  are  shown  in 
Appendix  A.  Distance  dx  and  d2  are  calculated  by  using  Equation  (20).  These  two 
distances  are  then  compared  by  calculating  absolute  differences  (d2-dx)  and  relative 
differences  {\^dl-d2\/dl).  However,  before  comparing  these  distances’  differences, 

results  of  distance  dx  and  d2  are  analyzed  statistically  as  shown  in  Table  2. 
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Phase  One:  Dive  #5 


Figure  14.  An  Example  Result  of  Phase  One:  Downward  Looking  View  of  DR 

and  DRNcom  Trajectories. 


Table  2.  Statistical  Analysis  of  Distance  dx  and  cl2 . 


Distance 

d \  (Pgps  "  P dr) 

d 2  (P gps  ~  P ncom) 

Mean 

168.895  m 

62.970  111 

Minimum 

38.421  m 

7.219  m 

Maximum 

295.417  m 

264.922  m 

Median 

167.170  m 

50.697  m 

Standard  Deviation 

66.147  m 

40.196  m 

Since  Ti  was  estimated  by  the  traditional  DR  algorithm  without  considering  the 
ocean  current  effects,  the  values  of  dx  were  larger  than  those  of  d2 .  Next,  the  absolute  and 

relative  improvements  of  To  compared  to  Ti  were  calculated  for  each  dive.  As  shown  in 
Figure  15,  by  applying  the  Regional  NCOM  prediction  data,  most  of  the  T2  trajectories 
improved  by  97.18%  compared  to  the  Ti  trajectories.  These  results  were  similar  to  those 
of  previous  work  from  Smith  et  al.  (2010a). 
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Absolute  Error:  d_  -  d„  Relative  Error:  (d„  -  dj  /  d, 

2  1  1  2'  1 


Figure  15.  Absolute  and  Relative  Differences  between  di  and  <i, . 

In  addition,  the  statistical  analysis  of  both  absolute  and  relative  differences 
between  dl  and  d2  was  provided  in  Table  3.  The  mean  of  absolute  difference  was 
-105.93  meters,  and  its  standard  deviation  was  69.21  meters.  The  root  mean  square  (RMS) 
of  distance  difference  was  126.25  meters.  Furthermore,  the  T2  improved  by  58.0%  in 
terms  of  mean  relative  distance  difference,  25.9%  in  terms  of  standard  deviation,  and 
63.4%  in  terms  of  RMS.  These  results  are  enough  to  account  for  the  ocean  current  effect 
in  estimating  a  glider’s  underwater  trajectory.  Based  on  these  results,  Phase  Two 
attempted  to  develop  a  novel  method  for  estimating  a  glider’s  optimal  underwater 
trajectory  by  applying  the  Regional  NCOM  to  the  DR  trajectory. 


Table  3.  Statistical  Analysis  of  the  Absolute  and  Relative  Differences 

between  di  and  dn . 


Differences 

Absolute  Values 

Relative  Ratio 

Improvement 

97.06  % 

Mean 

-105.93  m 

0.5799 

Standard  Deviation 

69.21  m 

0.2588 

RMS 

126.25  m 

0.6343 
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2. 


Results  of  Phase  Two 


This  section  discusses  the  results  of  applying  two  correction  methods  to  the  DR 
(Ti)  trajectory  in  four  steps.  In  the  first  two  steps,  two  different  corrected  trajectories,  the 
DR  trajectory  corrected  by  the  depth-averaged  correction  method  (C i)  and  the  DR 
trajectory  corrected  by  the  depth-dependent  correction  method  (C2),  were  estimated. 
Then,  to  compare  these  two  correction  methods,  an  optimal  trajectory  (C3)  was  estimated 
by  using  the  iteration  method  from  C2  until  the  relative  error  converges  to  less  than  10-'’ 
with  the  Regional  NCOM  prediction  data.  Finally,  the  Ci  and  C2  were  compared  to  C3  to 
determine  which  correction  method  is  better  for  estimating  a  glider’s  underwater 
trajectory  more  precisely. 

a.  Step  One:  DR  Trajectory  Corrected  by  the  Depth-Averaged  Correction 
Method 

Step  One  provides  the  results  of  the  depth-averaged  correction  method  when 
applied  to  the  DR  trajectory  (Ti).  Generally,  since  this  method  corrected  trajectory 
linearly,  the  shape  of  the  corrected  trajectory  (C 1)  looked  similar  to  T 1 .  Moreover,  the 
length  of  the  Ci  was  shorter  than  the  original  Ti.  Figure  16  shows  an  example  of  the 
results,  and  more  figures  of  results  are  shown  in  Appendix  A. 

b.  Step  Two:  DR  Trajectory  Corrected  by  the  Depth-Dependent  Correction 
Method 

Step  Two  provides  results  of  the  depth-dependent  correction  method  when 
applied  to  the  DR  trajectory  (T 1)  (see  Figure  17).  Actually,  this  step  was  supposed  to  deal 
with  the  DRNcom  trajectory  (T2),  but  the  results  of  this  correction  method  for  Ti  and  T2 
are  theoretically  the  same  because  T2  is  updated  from  Ti  by  considering  the  depth- 
dependent  Regional  NCOM  prediction  data.  Thus,  Ti  could  be  used  directly  without  any 
additional  processing  to  compensate  for  the  ocean  current  effects,  which  was  done  in 
Phase  One.  Therefore,  T2  did  not  need  to  be  corrected  with  either  of  the  two  correction 
methods  in  this  study.  The  mathematical  proof  is  presented  in  Appendix  B  due  to  its 
complexity. 
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Step  One:  Dive  #50 


Figure  16.  The  Result  of  Step  One:  Downward  Looking  View  of  the  Ti,  the  T3, 
and  the  Ci  with  the  Depth-averaged  Correction  Method. 


Step  Two:  Dive  #50 


Figure  17.  The  Result  of  Step  Two:  Downward  Looking  View  of  the  Ti,  the  T3, 
and  the  C2  with  the  Depth-dependent  Correction. 
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Unlike  the  depth-averaged  correction  method,  the  depth-dependent  correction 
method  corrected  trajectory  nonlinearly  using  the  Regional  NCOM  prediction  data, 
which  depends  on  the  depth,  time-consequent  location.  Since  this  method  corrected  Ti, 
the  shape  of  the  corrected  trajectory  (C2)  might  look  similar  to  Ti,  but  C2  was  more 
meandering  than  Ci  due  to  the  dependence  of  the  Regional  NCOM  prediction  data.  In 
addition  to  the  example  result  of  Step  Two  that  is  shown  in  Figure  17,  more  figures  of 
results  are  shown  in  Appendix  A. 

Both  Ci  and  C2  appear  to  be  more  accurate  than  the  original  Ti,  but  this  result  is 
inconclusive.  It  cannot  be  detennined  which  trajectory,  Ci  or  C2,  was  estimated  more 
precisely  due  to  the  lack  of  a  glider’s  actual  trajectory  information.  At  this  time,  results 
suggest  that  C2  was  estimated  better  than  T 1  using  the  depth-averaged  correction  method 
by  incorporating  the  Regional  NCOM  prediction  data.  It  is,  however,  difficult  to 
detennine  which  trajectory  was  more  precisely  compared  to  the  real  trajectory. 

c.  Step  Three:  Optimal  Trajectory  Development  by  Using  Iteration  Method 
from  the  C2 

To  address  the  problems  caused  in  the  previous  two  steps,  an  optimal  trajectory 
was  developed  as  a  reference  point  for  comparison  to  detennine  a  better  correction 
method.  It  is  not  an  easy  way  to  estimate  an  optimal  trajectory  because  there  is  less 
information:  the  linearly  interpolated  trajectory  (T3),  which  is  not  a  real  trajectory, 
between  each  GPS  fixed  position,  the  corrected  DR  trajectory  with  the  depth-dependent 
correction  method  (C2),  and  the  Regional  NCOM  prediction  data  as  the  ocean  currents 
information.  Furthermore,  since  the  Regional  NCOM  prediction  data  is  not  exactly  same 
with  the  real  ocean  currents  velocity,  it  is  assumed  that  the  ocean  currents  structure  of  the 
Regional  NCOM  is  similar  to  that  of  the  in-situ.  Given  these  data  and  assumption,  an 
optimal  trajectory  was  estimated  from  C2  by  using  the  iteration  method  (see  Figure  18). 
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Step  Three:  Dive  #50 


Figure  18.  The  Result  of  Step  Three:  Downward  Looking  View  of  C2  and  its 

Iterated  Trajectory,  C3. 

In  this  section,  the  relative  error  calculation  method  and  the  procedure  of  the 
iteration  method  will  be  explained  briefly.  The  details  of  the  iteration  method  are 
presented  in  Appendix  C,  and  more  figures  of  results  are  shown  in  Appendix  A. 

First,  the  relative  error  calculation  method  is  explained  as  follows.  Let  the 
reference  and  iterated  trajectory  be  Tm,  Tm+i,  the  sum  of  distance  difference  between  two 
trajectories  be  Dm ,  and  the  total  length  of  the  reference  trajectory  be  Lm  at  each  step. 

Then  the  relative  error  ( Erel )  at  each  step  can  be  calculated  by  dividing  Dm  by  Lm  as 
shown  in  Equation  (33). 

Erel  =  Dm/Lm,  m  =  l,  -\,M  (33) 

Now,  the  procedure  of  the  iteration  method  will  be  explained.  In  the  first  step,  it  is 
assumed  that  the  initial  reference  trajectory  is  T3,  and  then  the  C2  is  estimated  from  Ti  by 
applying  the  depth-dependent  correction  method  against  the  distance  difference  between 
P dr  and  PgpS  ( AXW ,  see  Figure  4).  The  used  Regional  NCOM  prediction  data  are  obtained 
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from  each  point  on  T3.  After  estimating  the  C2,  the  relative  error  is  computed  by  using 
Equation  (33),  but  this  step  is  not  included  in  the  iteration  step. 

After  the  first  step,  C2  becomes  a  reference  trajectory,  and  the  Regional  NCOM 
prediction  data  is  obtained  from  each  point  on  C2  as  in  the  first  step.  Then  Ti  is  corrected 
again  by  using  the  depth-dependent  correction  method  against  the  same  distance 
difference  AT11  as  the  first  step.  This  time  is  regarded  as  the  first  iteration.  After  the  first 
iteration,  the  relative  error  is  calculated.  If  the  value  of  relative  error  is  less  than  1 0-\  the 
iteration  is  finished,  and  the  final  trajectory,  which  is  estimated  by  the  iteration  method, 
would  be  the  C3. 

In  this  study,  the  mean  iteration  steps  were  3.088  with  the  maximum  of  eight  steps, 
and  the  mean  relative  error  at  Step  One,  Two,  and  Three  was  0.387,  0.001,  and  0.003, 
respectively.  Table  4  provides  the  statistical  results  of  the  iteration  method. 


Table  4.  Statistical  Analysis  of  the  Iteration  Method. 


Relative  Error 

First  Step 

Second  Step 

Third  Step 

Mean 

0.387 

0.001 

0.003 

Minimum 

0.002 

1.309xl0'6 

0.000 

Maximum 

20.056 

0.448 

0.193 

Median 

0.034 

0.059 

0.023 

Standard  Deviation 

2.431 

9.646*  10'5 

3.972xl0'7 

Even  though  an  optimal  trajectory  can  be  estimated  with  iteration  method,  the 
precision  of  the  trajectory  is  not  guaranteed.  For  the  purpose  of  this  study,  however,  it  is 
assumed  that  C3  is  the  most  similar  trajectory  to  a  real  trajectory. 

d.  Step  Four:  Comparison  of  Three  Types  of  Corrected  Trajectories 

This  section  provides  results  of  comparison  between  two  types  of  corrected 
trajectory,  Ci  and  C2,  to  an  estimated  optimal  trajectory,  C3,  to  detennine  which 
correction  method  is  better  for  estimating  a  glider’s  optimal  underwater  trajectory.  Figure 
19  presents  comparison  of  these  three  trajectories,  and  more  figures  of  results  can  be  seen 
in  Appendix  A. 
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Step  Four:  Dive  #50 


Figure  19.  Downward  Looking  View  of  Three  Types  of  Corrected  Trajectory. 

Given  these  three  trajectories,  the  RMS  of  distance  difference  between  Ci  and  C3 
( d,  ),  C2  and  C3  (d4)  for  each  dive  was  calculated  to  determine  a  better  method.  Table  5 

provides  statistical  analysis  of  the  RMS  of  distance  d,  [RMS  ( d3 )]  and  d4  [RMS  ( dA )] 
for  total  dives,  and  Figure  20  shows  absolute  and  relative  differences  between  RMS  ( d3  ) 
and  RMS  ( d4  ). 


Table  5.  Statistical  Analysis  of  the  RMS  ( d3  )  and  RMS  ( d4  )  for  Total 

Dives. 


RMS  of  Distance 

d3  (Ci  -  C3) 

d4  (C2  -  C3) 

Mean 

9.890  m 

1.034  m 

Minimum 

1.195  m 

0.003  m 

Maximum 

80.536  m 

59.100  m 

Median 

6.316  m 

0.052  m 

Standard  Deviation 

11.052  m 

7.168  m 
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From  the  Table  5,  all  values  of  RMS  (d4)  are  less  than  RMS  ( d,  ),  which  means 

the  C2  was  more  close  to  the  C3.  In  addition,  from  Figure  20,  it  can  be  said  that  the  C2 
was  better  estimated  than  Ci  by  comparing  the  absolute  and  relative  differences  between 
RMS  (d3)  and  RMS  ( d4 ). 


Absolute  Differences:  RMS(d4)  -  RMS(d3)  for  Total  Dives 


Relative  Differences:  [RMS(d3)  -  RMS(d4)]  /  RMS(d3)  for  Total  Dives 


Figure  20.  Absolute  and  Relative  Differences  between  RMS  ( d3  )  and 

RMS  ( d4  )  for  Total  Dives. 


Table  6  provides  statistical  analysis  of  the  maximum  distance  of  d3  and  d4  for 
total  dives.  Like  the  analysis  of  RMS  of  distance,  all  values  of  d3  are  larger  than  those  of 
d4  .  This  also  means  C2  was  closer  to  C3. 


Table  6.  Statistical  Analysis  of  the  Maximum  Distance  of  d  ,  and  d4  for 

Total  Dives. 


Maximum  Distance 

d3  (Ci  -  C3) 

d4  (C2  -  C3) 

Mean 

16.052  m 

2.142  m 

Minimum 

1.912  m 

0.005  m 

Maximum 

142.648  m 

123.491  m 

Median 

10.070  m 

0.095  m 

Standard  Deviation 

19.354  m 

14.987  m 
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In  addition  to  the  statistical  analysis,  Figure  21  describes  the  maximum  distance 
of  d3  and  d4  for  total  dives.  Only  four  of  the  data  of  d4  were  longer  than  one  meter,  but 

no  data  of  d3  was  less  than  one  meter. 


Maximum  Distance  of  d,  and  d,  for  Total  Dives 

3  4 


Figure  21.  Maximum  Distance  of  d3  and  d4  for  Total  Dives. 

By  comparing  these  results,  the  tentative  conclusion  can  be  expressed  that  C2, 
which  was  corrected  by  the  depth-dependent  correction  method,  was  estimated  more 
precisely  than  Ci,  which  was  corrected  by  the  depth- averaged  correction  method. 

Now,  a  controversial  issue  can  be  argued  as  to  which  trajectory  was  estimated 
more  precisely  between  C2  and  C3.  It  is  expected  that  C3  would  be  more  precise  than  C2 
because  C3  was  estimated  by  the  iteration  method  by  considering  the  Regional  NCOM 
prediction  data,  which  depends  on  the  depth  and  time-consequent  location,  on  the 
corrected  trajectories.  Even  though  the  iteration  method  makes  the  trajectory  more 
precise  than  without  iteration,  the  depth-dependent  correction  method  would  be  more 
efficient  than  the  iteration  method  because  the  iteration  method  is  complex  and  difficult 
to  deal  with,  so  it  would  take  more  processing  time  to  estimate  C3  requires  than  that  to 
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estimate  C2.  Furthermore,  the  distance  difference  between  C2  and  C3  was  minimal. 
Therefore,  if  the  time  cost  of  data  processing  is  larger  than  the  distance  difference  cost,  it 
is  expected  that  the  depth-dependent  correction  method  without  combining  with  the 
iteration  method  would  be  an  efficient  and  precise  enough  method  for  estimating  a 
glider’s  optimal  underwater  trajectory. 

B.  DISCUSSION 

This  section  discusses  the  relevance  of  findings  in  this  study  to  the  ROKN.  It  is 
divided  into  three  parts:  advantages  of  an  underwater  glider  in  ASW,  relevance  of 
findings  to  ASW,  and  relevance  to  future  ASW  in  the  ROKN.  The  first  section  describes 
three  high-level  advantages  of  the  underwater  glider  in  ASW:  cost-effectiveness,  various 
mission  capabilities,  and  risk  reduction.  The  second  section  discusses  relevance  of 
findings  to  ASW  by  reviewing  the  U.S.  Navy  Unmanned  Undersea  Vehicle  Master  Plan 
(UUVMP)  (DON  2004)  and  by  applying  four  Sub-Pillar  missions  to  the  underwater 
glider  in  ASW.  An  optimally  estimated  trajectory  would  be  beneficial  to  operating  ISR  in 
the  area  of  interest,  detecting  and  identifying  the  enemy’s  bottoming  submarine, 
enhancing  and  supporting  tactical  oceanography  by  providing  spatiotemporally  high 
resolution  ocean  data,  and  providing  communication  and  localization  infonnation. 
Finally,  we  discuss  the  relevance  of  these  findings  to  future  ASW  by  reviewing  ASW 
threats  to  the  ROKN  and  by  suggesting  recommendations  for  using  the  underwater  glider 
for  ASW  in  the  ROKN. 

1.  Advantages  of  the  Underwater  Glider  in  ASW 

An  underwater  glider  has  primarily  been  used  for  oceanographic  research  with 
data  collection  or  environment  monitoring.  Recently,  however,  there  is  an  increasing 
tendency  to  take  advantage  of  the  glider  in  naval  operations,  especially,  ASW.  Some 
UUVs  were  already  used  in  previous  naval  operations;  for  example,  the  REMUS  was 
employed  for  mine  countenneasures  (MCM)  operations  during  Operation  Iraqi  Freedom 
(Clegg  and  Peterson  2003).  However,  the  glider  has  not  been  deployed  in  real-time 
combat  operations  yet,  and  the  development  of  the  glider’s  capability  and  concept  of 
operations  (CONOPs)  for  ASW  is  still  ongoing.  Even  though  the  glider  needs  to  be 
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improved  in  various  fields,  it  is  expected  that  the  glider  will  be  an  effective  platform  in 
the  future  ASW  due  to  their  various  high-level  advantages  as  follows. 

a.  Cost  Effectiveness 

An  underwater  glider  is  designed  for  long-duration  and  long-range  missions.  Its 
propulsion  system  requires  less  energy,  and  it  is  cheaper  than  ships  or  other  marine 
vehicles,  which  perform  similar  missions.  Thus,  its  low  cost  allows  multiple  gliders  to  be 
deployed  and  to  patrol  in  a  broader  area  than  other  ships.  Furthennore,  they  can  collect 
more  spatiotemporally  vast  ocean  data  than  existing  ocean  observation  ships  or 
instruments.  For  these  reasons,  the  underwater  glider  is  cost-effective  platform. 

b.  Various  Mission  Capabilities 

A  glider  can  be  equipped  with  additional  sensors  for  various  missions;  for 
example,  acoustic  sonar  for  ISR,  side  scan  sonar  (SSS)  for  detecting  a  bottoming 
submarine,  and  some  optic  sensors  for  detecting  and  identifying  unknown  underwater 
targets.  Even  though  installing  additional  sensors  shortens  a  glider’s  mission  time,  the 
navy  could  benefit  from  a  glider’s  mission  capabilities.  This  is  because  a  glider  can  cover 
wider  mission  area  and  perform  for  a  longer  time  with  less  cost  than  presently  used  ships 
or  other  marine  vehicles.  In  addition,  the  glider’s  cost-effectiveness  allows  multiple 
gliders  to  be  deployed  simultaneously.  For  example,  a  fleet  of  gliders  would  be  able  to 
act  as  an  intelligence  network  by  gathering  underwater  target  information  simultaneously. 
By  doing  so,  they  can  provide  more  precise  target  infonnation,  so  that  the  ASW 
capability  would  be  increased.  Moreover,  multiple  gliders  can  act  as  communication 
nodes  with  other  platfonns  such  as  surface  ships,  submarines,  and  other  UUVs  in  the 
subsurface  area.  For  these  reasons,  a  glider’s  ASW  capability  effectiveness  would  be 
increased. 


c.  Risk  Reduction 

Even  though  an  underwater  glider  has  various  mission  capabilities,  it  is  still 
required  to  deploy  manned  platfonns  to  accomplish  the  ROKN’s  mission  more 
effectively.  Manned  platforms,  however,  have  an  ever-present  danger  when  they  perfonn 
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missions  in  a  hazardous  operation  environment  like  a  shallow-water,  mine  burial  area. 
Since  an  underwater  glider  is  an  unmanned  vehicle  and  is  operated  autonomously, 
manned  platfonns  can  be  replaced  by  the  glider  in  such  perilous  AO;  for  example, 
explosive  ordnance  disposal  (EOD)  and  mine  localization.  Especially,  the  glider  could  be 
a  safe  platform  when  perfonning  ASW  operations  because  it  is  relatively  quieter  than 
other  UUVs  due  to  its  buoyancy-driven  engine.  The  quietness  makes  the  glider  generate 
less  self-noise  so  it  has  superior  detection  ability  when  it  attempts  to  detect  underwater 
targets  by  using  acoustics.  Furthennore,  the  quietness  makes  the  enemy  have  more 
difficulty  detecting  a  glider.  Therefore,  an  underwater  glider  would  be  able  to  perfonn  its 
mission  more  safely  and  it  offers  more  opportunity  to  accomplish  its  mission  more 
successfully. 

2.  Relevance  of  Findings  to  ASW 

The  following  two  subsections  discuss  the  relevance  of  findings  from  this  study  to 
ASW  in  the  ROKN.  The  first  subsection  outlines  the  overview  of  the  U.S.  Navy  UUVMP 
to  provide  the  background  of  the  UUV’s  nine  Sub-Pillar  missions  based  on  the  Sea  Power 
21  pillars.  The  second  subsection  discusses  the  relevance  of  findings  to  ASW  related  to 
the  UUVs’  mission  capabilities  mentioned  in  the  Master  Plan.  This  section  constrains  the 
area  of  cover  to  four  capabilities  related  to  ASW:  ISR;  Littoral  ASW;  Tactical 
Oceanography;  and  Communication/Navigation  Network  Nodes  (CN3). 

a.  Overview  of  the  U.S.  Navy  UUVMP 

The  U.S.  Navy  UUVMP  (hereafter  referred  to  as  UUVMP)  is  an  updated 
UUVMP  based  on  the  2000  UUVMP.  It  was  updated  due  to  changes  in  various  fields  in 
the  U.S.  Navy:  technology,  platforms,  and  other  factors  (DON  2004).  The  UUVMP 
classifies  the  high-priority  missions  based  on  the  four  Sea  Power  2 1  pillars:  Sea  Shield, 
Sea  Strike,  Sea  Base,  and  FORCEnet  (DON  2004).  Then,  these  four  pillars  are  divided  by 
UUV  missions  into  nine  sub-pillar  missions  as  shown  in  Figure  22,  and  the  details  of 
each  mission  are  described  in  the  Master  Plan.  However,  not  all  sub-pillar  missions  can 
be  performed  by  an  underwater  glider  due  to  its  minimal  velocity  and  limited  payload. 
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For  these  reasons,  we  will  discuss  the  relevance  of  findings  related  to  only  some  of  those 
capabilities  in  ASW  operations:  ISR;  Littoral  ASW;  Tactical  Oceanography;  and  CN3. 


Figure  22.  Sea  Power  21  Nine  Sub-Pillar  Capabilities  (from  DON  2004). 


b.  Relevance  Related  to  the  UUVMP 

In  this  section,  the  relevance  of  findings  to  the  ROKN  that  relate  to  the  UUVMP 
and  ASW  operations  will  be  detailed  in  four  sub-pillar  missions  as  follows. 

(1)  ISR 

ISR  is  an  essential  operation  in  ASW.  The  ISR  mission  consists  of  intelligence 
gathering,  target  detection  and  localization,  and  mapping  (DON  2004).  An  underwater 
glider  can  be  used  for  persistent  and  clandestine  ISR,  and  it  can  perform  the  mission 
effectively  with  high-endurance  and  long-range  capabilities,  its  ability  to  operate  in 
shallow  water  and  coastal  areas,  and  its  autonomous  features  (DON  2004).  The  glider  can 
patrol  a  friendly  harbor  or  coastal  area  near  the  high-value  infrastructures.  Even  though 
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the  glider  is  limited  in  its  ability  to  attack  an  enemy  submarine  directly  due  to  its  minimal 
velocity  and  limited  payload,  it  can  provide  more  precise  and  reliable  target  information 
in  the  area  of  interest.  Thus,  the  ROKN  can  get  the  target  information  more  precisely  and 
prevent  the  enemy  from  infiltrating  into  friendly  territory.  This  precise  information  makes 
the  navy  spend  less  effort  to  detect,  to  identify,  and  to  attack  the  enemy  underwater  forces 
in  ASW. 

(2)  Littoral  ASW 

The  U.S.  Navy  Task  Force  ASW  has  instituted  a  new  focus  on  littoral  ASW  and 
described  three  major  categories  in  ASW  as  follows  (DON  2004): 

“Hold  at  Risk”:  monitoring  all  the  submarines  that  exit  a  port  or  transit  a 
checkpoint,  “Maritime  Shield”:  clearing  and  maintaining  a  large  Carrier  or  Expeditionary 
Strike  Group  (CSG  or  ESG)  operation  area  free  of  threat  submarines,  and  “Protected 
Passage”:  clearing  and  maintaining  a  route  for  an  ESG  from  one  operating  area  to  another 
free  of  threat  submarines. 


Figure  23.  Schema  of  Task  Force  ASW  Three  Major  Categories  (from  DON 

2004). 

An  underwater  glider  can  offer  some  support  in  the  second  and  third  categories 
but  is  more  suitable  for  “Hold  at  Risk”  missions  due  to  the  glider’s  restricted  velocity 
(DON  2004).  The  glider  can  perform  the  “Hold  at  Risk”  mission  more  effectively  by 
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connecting  communication  with  other  platforms  and  by  providing  longer  detection  range 
(DON  2004),  estimating  its  location  precisely  consequent  the  enemy  location.  Moreover, 
since  the  glider  can  operate  in  shallower  water  and  closer  to  shore  than  other  platfonns,  it 
can  establish  an  early  warning  barrier.  This  early  warning  barrier  allows  friendly  manned 
platforms  to  stay  in  a  safe  area  away  from  the  enemy’s  detection  or  hazardous 
environment  and  to  focus  on  preparing  the  next  mission  after  the  glider  detects  and 
identifies  some  underwater  objects.  Even  though  the  UUVMP  analyzed  and  determined 
to  be  suitable  for  this  “Hold  at  Risk”  mission  was  a  large  vehicle  class  of  UUV,  we  think 
that  an  underwater  glider  also  would  be  able  to  accomplish  the  mission  with  some 
limitations. 

In  addition  to  the  “Hold  at  Risk”  mission,  a  glider  would  be  able  to  detect  and 
identify  the  enemy  underwater  forces  when  they  are  sitting  on  the  seabed.  Such  a  mission 
could  be  performed  more  effectively  by  using  SSS-like  mine  localization.  The  detected 
target  infonnation  would  be  more  precise  when  a  glider’s  position  consequent  the  target 
position  is  estimated  precisely.  By  doing  so,  a  commander  can  evaluate  the  target  quickly 
and  make  the  right  decision  so  that  ASW  operation  would  be  performed  effectively. 

(3)  Tactical  Oceanography 

Knowing  the  environment  feature  of  AO  is  an  essential  component  in  ASW.  An 
underwater  glider  can  collect  spatiotemporally  high-resolution  ocean  data  during  a  long¬ 
term  mission,  and  these  data  provide  maritime  domain  awareness  (MDA)  (DON  2007) 
and  intelligence  preparation  of  the  operational  environment  (IPOE).  These  collected  data 
can  be  used  for  increasing  ocean  model  prediction  precision  and  antisubmarine  warfare 
tactical  decision  aids  (TDAs)  to  support  and  to  optimize  operational  planning  and  asset 
management  (DON  2004).  Thus,  the  location  where  these  data  are  collected  is  very 
important.  These  spatiotemporally  high-resolution  ocean  data  can  be  obtained  by 
calculating  from  the  glider’s  precise  underwater  trajectory.  For  example,  CTD  and  some 
optical  data,  which  are  collected  accurately  and  persistently,  would  optimize  the 
oceanographic  characteristic  of  AO  to  help  with  understanding  environmental  conditions 
that  impact  the  operational  effectiveness  in  ASW. 
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(4)  CN3 

An  underwater  glider  has  restricted  CN3  capabilities  due  to  its  payload  limitation. 
Operators,  however,  still  can  take  some  advantage  of  the  glider’s  CN3  capabilities. 
Communication  is  an  essential  factor  in  battle,  but  it  is  limited  in  a  subsurface  mission 
area  due  to  the  acoustic  wave’s  characteristics:  attenuation,  refraction,  and  reflection. 
Although  a  glider  is  limited  in  providing  long-range  communication  between  subsurface 
forces,  it  can  act  as  a  network  node  as  previously  mentioned  in  Section  1.  Furthermore, 
multiple  gliders  can  act  as  communication  nodes  with  other  platforms  such  as  surface 
ships,  submarines,  and  other  UUVs  in  the  subsurface  area. 

An  even  more  important  feature  is  that  an  underwater  glider  can  provide 
underwater  target  location  infonnation.  This  target  location  infonnation  would  be  precise 
when  the  glider’s  underwater  trajectory  is  estimated  precisely.  In  addition,  a  fleet  of 
gliders  would  be  able  to  gather  underwater  target  infonnation  simultaneously.  Thus,  they 
can  provide  more  precise  target  infonnation  than  when  a  single  glider  operates.  For  these 
reasons,  a  glider’s  ASW  capability  effectiveness  would  be  increased. 

3.  Relevance  to  Future  ASW  in  the  ROKN 

This  section  discusses  relevance  to  future  ASW  in  the  ROKN  by  reviewing  ASW 
threats  to  the  ROKN  and  providing  future  recommendations  for  ASW  in  the  ROKN.  The 
first  section  describes  the  force  power  and  organization  of  the  North  Korea  Navy’s 
underwater  forces,  which  are  the  most  threatening  underwater  forces  to  the  ROKN.  The 
second  section  provides  some  future  recommendations  for  ASW  in  the  ROKN  to 
overcome  present  limitations  in  the  ASW  platfonn  aspect. 

a.  ASW  Threats  to  the  ROKN 

Even  though  ASW  operation  was  conducted  by  the  ROKN  a  bit  later  than  by 
other  navies,  it  has  been  become  an  important  naval  operation  since  the  Korean  War.  The 
most  threatening  forces  to  the  ROKN  are  the  North  Korea  Navy’s  underwater  forces.  The 
North  Korea  Navy  has  about  70  underwater  forces,  which  consist  of  R  class  and  Sang-0 
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class  submarines  and  Yono  class  submersible  vessels,  as  shown  in  Figure  24 
(Government  of  the  Republic  of  Korea  2011;  MND  2013). 


(a)  (b)  (c) 


Figure  24.  Underwater  Forces  of  The  North  Korea  Navy:  (a)  R  class 

Submarine;  (b)  Sang-0  class  Submarine;  (c)  Yono  class  submersible 
vehicles  (from  Government  of  the  Republic  of  Korea,  2011). 

Their  mission  is  a  sea  line  of  communication  (SLOC)  interdiction,  mine  laying, 
attack  surface  ships,  and  supporting  the  infiltration  of  special  warfare  units  (MND  2013). 
In  fact,  North  Korea  has  been  continuing  several  underwater  infiltrations  and 
provocations  since  the  Korean  War.  Recently,  on  March  26,  2010,  in  a  tragic  incident,  the 
ROK  Ship  Cheonan  was  sunk  in  a  surprise  torpedo  attack  by  North  Korea  Navy’s  small 
submarine  (MND  2010).  In  addition,  North  Korea  continues  making  and  developing  new 
submarines  and  underwater  weapons  in  order  to  augment  the  underwater  forces  as 
asymmetrical  capabilities  (MND  2013).  Therefore,  such  underwater  forces  of  North 
Korea  pose  the  greatest  threat  to  the  ROKN’s  ASW  operation. 

The  lessons  from  the  review  of  ASW  threat,  the  potential  power  of  submarines, 
and  the  importance  of  effort  and  resources  in  ASW  to  counter  them  were  learned.  ASW 
will  not  disappear  as  long  as  enemy  submarines  are  present,  and  they  will  continue  to  be 
developed  and  pose  the  most  threats.  Therefore,  the  ROKN  who  is  faced  with  these 
threats  must  develop  active  and  forward-looking  countermeasures  in  order  not  to  repeat 
the  failures  which  were  demonstrated  in  the  past. 

b.  Future  Recommendations  to  ASW  in  the  ROKN 

An  effective  platform  development  could  be  considered  as  a  part  of  improving 
ASW  capability  in  the  ROKN.  Currently,  various  types  of  ASW  platforms  are  being 
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operated  in  the  ROKN:  surface  ships,  submarines,  aircrafts  (see  Figure  25).  These 
platforms,  however,  have  some  limitations — such  as  shortage  of  platform  numbers,  time 
available  of  platform,  and  limited  manpower  and  budget  compared  to  AO — that  must  be 
overcome.  Since  an  underwater  glider  has  various  advantages  as  mentioned  in  Section  1, 
it  can  be  used  to  compensate  for  these  limitations.  Thus,  an  underwater  glider  is  expected 
to  be  applied  to  conduct  effective  ASW  operations  in  the  near  future.  The  types  of  glider 
could  be  selected  depending  on  detailed  missions  and  the  environment.  If  a  glider  is  to  be 
used  for  long-term  missions  like  collecting  ocean  data  and  conducting  ISR  near  the  deep 
and  open  ocean,  the  Spray  or  the  Seaglider,  which  are  optimized  for  deep-water 
operations,  would  be  a  proper  platform.  On  the  other  hand,  if  a  glider  is  to  be  used  for 


short-term  missions  like  littoral  ASW,  including  ISR  near  the  shallow  or  coastal  ocean, 
the  Slocum  Electric,  which  is  optimized  for  shallow-water  operations,  would  be  a  suitable 


platfonn. 


Figure  25.  Currently  Operating  Platforms  for  ASW  in  the  ROKN:  (a)  Surface 
Ship  (DDG)  (from  MND  2014,  https://www.flickr.com/photos 
/konnnd/74459694 1 6/in/photolist-fq  1 4Rx-fq  1 4Hp-fq  1 4T4-fq  1 5 1 H- 
ckWcsd-NAzrw-ckYxi3-ktTP6w-9sTYsu-ajkxJd-mB5CqH-cHudGY- 
dU98Hp-dUeKiU-dUeKpE-dU98yV-dUeKt5-dU98De-dUeKkm- 
498nUX);  (b)  Submarine  (214  Class)  (from  Center  for  International 
Maritime  Security  [CIMSEC]  2014,  http://cimsec.org/wp-content/ 
uploads/20 14/02/Type_2 14_Conventional_Attack_Submarine- 
672x372.jpg);  (c)  Antisubmarine  Patrol  Aircraft  (P-3C)  (from  Aviation 
WA  2014,  http://www.aviationwa.org.au/wp-content/uploads 
/20 14/05/20 1404 1 7_950905_Lockheed_P-3C- 
III+_Orion_Keith_Anderson_2.jpg);  (d)  Antisubmarine  Helicopter 
(LYNX)  (from  the  ROKN  2014,  http://cfile25.uf.tistory.com/image 
/177CF81C4B1F334A2B8922). 


Future  ASW  CONOPs  would  be  network  centric  warfare  (NCW)  based  on 
Command,  Control,  Communications,  Computers,  and  ISR  (C4ISR).  Even  though  a 
present  glider  is  equipped  with  telecommunication  systems,  its  capability  is  limited 
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during  the  subsurface  time.  However,  if  a  fleet  of  gliders  fonns  an  intelligence  network, 
they  would  be  able  to  provide  more  precise  and  useful  infonnation.  Additionally,  if  there 
is  no  network  sensor  near  a  mission  area,  a  glider  can  act  as  a  network  node  itself 
between  the  glider  and  both  the  subsurface  and  surface  platforms.  Consequently,  ASW 
operation  capability  would  be  doubled  with  the  development  of  the  communication 
system  between  a  glider  and  operators  and  other  types  of  platforms. 

Such  a  beneficial  underwater  glider  would  be  an  effective  future  ASW  platform  in 
the  ROKN.  The  ROKN  is,  however,  in  the  early  stages  of  developing  its  glider. 
Moreover,  operating  the  glider  is  not  simple  compared  to  other  ocean  observation 
instruments.  Unlike  these  other  instruments,  the  glider  considers  the  operation 
environment  and  requires  more  effort  from  its  operators  to  receive  collected  data  and  to 
command  missions.  Thus,  the  ROKN  needs  to  develop  CONOPs,  as  well  as  manpower 
for  operation  and  maintenance.  Therefore,  to  take  advantages  of  the  glider  in  ASW,  the 
ROKN  should  train  its  operators,  develop  its  infrastructure,  and  systematically  prepare 
for  operating  the  underwater  glider. 

C.  CHAPTER  SUMMARY 

In  this  chapter,  we  discussed  the  results  of  corrected  trajectory  by  applying  our 
two  correction  methods,  the  depth-averaged  and  the  depth-dependent  correction  method, 
and  the  iteration  method.  By  using  these  methods,  the  glider’s  underwater  trajectories 
were  corrected  and  its  optimal  trajectories,  which  act  as  reference  trajectories,  were 
obtained  as  well.  Then,  we  could  compare  corrected  trajectories  to  an  optimal  trajectory 
at  each  dive  to  detennine  which  correction  method  is  better.  From  the  comparison  among 
these  three  trajectories,  we  could  get  a  tentative  conclusion  that  the  novel  method,  the 
depth-dependent  correction  method,  is  more  efficient  and  precise  enough  to  estimate  the 
glider’s  optimal  underwater  trajectory. 

Based  on  these  finding  from  this  study,  we  discussed  the  relevance  of  findings  to 
the  ROKN  by  describing  advantages  of  the  underwater  glider  and  by  reviewing  the 
UUVMP.  The  underwater  glider  has  three  high-level  advantages  in  ASW:  cost 
effectiveness,  various  mission  capabilities,  and  risk  reduction.  With  these  advantages,  the 
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glider  can  accomplish  various  naval  operations,  especially  ASW  operations.  To  discuss 
the  relevance  of  these  findings  to  the  ROKN,  the  UUVMP  was  reviewed  to  get  the 
background  on  UUVs’  sub-pillar  capabilities.  However,  since  a  glider  is  limited  in 
velocity  and  payload,  only  four  capabilities,  i.e.,  ISR,  Littoral  ASW,  Tactical 
Oceanography,  and  CN3,  were  covered  in  this  study.  Even  though  a  glider  still  has  some 
limitations  in  performing  these  four  missions,  it  is  expected  that  a  glider  would  be  able  to 
contribute  to  increased  ASW  capability  effectiveness.  Finally,  we  reviewed  ASW  threats 
to  the  ROKN  to  propose  future  recommendations.  An  underwater  glider  would  be  an 
effective  platform  that  can  solve  the  limited  ASW  platfonn  problem  in  the  future  ROKN. 
A  glider’s  telecommunication  system  would  allow  the  glider  to  perform  more  effectively 
in  ASW  related  to  NOW.  A  glider  could  be  a  network  node,  and  a  fleet  of  gliders  would 
be  able  to  provide  more  precise  intelligence.  To  improve  efficiency  of  a  glider,  different 
types  of  gliders  must  be  considered  with  their  mission  type  and  the  ocean  environment.  In 
addition,  since  a  glider  is  not  developed  well  enough  to  apply  to  naval  operation  in  the 
ROKN,  professional  and  systematic  preparation  is  required  to  take  advantage  of  the 
underwater  glider  in  ASW. 
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V.  CONCLUSIONS  AND  FUTURE  RESEARCH 


This  chapter  concludes  the  thesis  by  summarizing  the  results  of  the  method 
developed  by  the  thesis,  by  discussing  limitations  of  this  study,  and  by  suggesting  future 
research. 

A.  CONCLUSIONS 

This  study  was  intended  to  develop  a  novel  method  for  estimating  a  glider’s 
optimal  underwater  trajectory,  considering  the  ocean  current  effects.  A  precisely 
estimated  trajectory  would  be  beneficial  for  ocean  observation  and  naval  applications. 
Estimating  a  glider’s  underwater  trajectory  is,  however,  a  challenging  task  because  a 
glider’s  navigational  accuracy  is  limited  during  the  subsurface  missions.  Its  DR  algorithm 
does  not  consider  the  ocean  current  effects,  and  the  submerged  glider  is  limited  in  its 
ability  to  receive  GPS  fix  infonnation.  In  addition,  a  glider  lacks  inexpensive  and 
efficient  onboard  sensors  to  measure  ocean  currents.  For  these  reasons,  developing  a  new 
cost-efficient  and  reliable  method  is  needed. 

Two  methods  were  developed  to  solve  this  challenging  problem  by  applying  the 
Regional  NCOM  data  to  the  traditional  DR  trajectory.  Two  corrected  trajectories  (Ci  and 
C2)  were  compared  to  an  optimal  trajectory  (C3),  which  was  estimated  by  using  the 
depth-dependent  correction  method  and  incorporating  the  iteration  method.  Since  no  one 
knows  the  real  trajectory,  it  was  assumed  that  C3  is  the  most  precisely  estimated 
trajectory.  This  trajectory,  however,  requires  more  time  than  C2  to  process  the  data  of  the 
glider  and  the  ocean  model  data,  and  the  iteration  method  itself  is  lengthy  and  difficult  to 
understand.  Moreover,  the  distance  difference  between  C2  and  C3  is  minimal.  Thus,  the 
depth-dependent  correction  method  without  using  the  iteration  method  would  be  a  more 
efficient  and  sufficiently  precise  correction  method. 

Given  that  the  depth-dependent  correction  method  is  reliable  in  estimating  a 
glider’s  optimal  underwater  trajectory,  an  underwater  glider  could  be  more  beneficial  for 
ocean  observation  and  naval  applications.  A  glider  can  collect  spatiotemporally  higher- 
resolution  ocean  data  by  performing  persistent  ocean  sampling,  and  the  collected  data 
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would  be  more  valuable  if  the  glider’s  estimated  underwater  positions  are  precise. 
Furthermore,  the  precisely  estimated  underwater  trajectory  of  a  glider  would  have  a  great 
impact  on  ASW;  for  example,  perfonning  ISR,  detecting  and  identifying  the  enemy  target, 
using  multiple  gliders  to  form  a  communication  system  by  acting  as  network  nodes  and 
providing  precise  intelligence,  and  supporting  tactical  oceanography. 

Although  a  glider  has  many  advantages  in  ASW,  operating  it  is  more  difficult  and 
complex  than  operating  other  ocean  monitoring  instruments.  Thus,  more  professional  and 
systematic  preparation  is  required  to  apply  a  glider  to  future  naval  operations. 

B.  LIMITATIONS  AND  FUTURE  RESEARCH 

To  further  the  research  presented  in  this  paper,  several  future  research  areas  are 
suggested.  First  of  all,  since  the  C3  was  not  an  actual  trajectory,  the  result  of  this  study 
was  a  tentative  conclusion.  Thus,  to  estimate  a  more  precise  trajectory  and  to  verify  our 
novel  method,  a  glider’s  real  underwater  trajectory  information  must  be  obtained  by 
using  additional  sensors  like  ADCP  or  DVL  (see  Figure  26).  Then,  we  can  compare  the 
C3  to  the  real  trajectory  to  determine  which  correction  method  is  more  suitable. 


(a)  (b) 


Figure  26.  (a)  Sentinel  ADCP  (from  Teledyne  RD  Instruments  2014, 

http://www.rdinstruments.com/images/web_sentinell  105.jpg);  (b) 

Explorer  DVL  (from  Teledyne  RD  Instruments  2014, 
http://www.rdinstruments.com/images/explorer_pa_pd2.jpg). 

Second,  this  study  used  limited  data.  The  Slocum  conducted  a  total  of  68  dives  for 
a  three-day  deployment,  and  the  length  of  each  dive  was  short,  with  362  meters  of 
average  straight  DR  trajectory  distance.  It  was  a  relatively  short  experiment  compared  to 
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previous  works.  For  example,  Smith  et  al.  (2010a)  used  the  Slocum  data,  which  were 
obtained  from  a  month-long  deployment  with  an  average  trajectory  length  of  two 
kilometers  and  more  than  200  times  of  surfacing.  However,  since  our  particular 
experiment  was  not  originally  intended  for  this  study,  it  was  uncontrollable.  Generally,  an 
underwater  glider  is  designed  for  long-range  and  long-duration  missions.  Thus,  we  need 
to  research  more  realistically  to  detennine  the  suitability  of  our  novel  method  by  using 
data  from  long-range  and  long-duration  missions. 

Third,  this  study  used  relatively  low  resolution  of  the  Regional  NCOM  data.  Our 
average  dive  length  was  362  meters,  but  the  resolution  of  the  Regional  NCOM  data  set 
was  1/30  degree  (about  3.6  kilometers).  Generally,  we  can  say  that  1/30  degree  is  high 
resolution  because  the  Global  HYCOM  has  1/12  degree  (about  nine  kilometers),  but  it 
was  actually  a  relatively  low  resolution  to  the  glider’s  dive  length.  Velocity  change  of  the 
ocean  currents  in  one  grid  cell  of  the  Regional  NCOM  data  was  not  sensitive,  and  each 
dive  length  was  included  in  one  grid  cell  of  the  Regional  NCOM  data.  Thus,  the  velocity 
effect  on  estimating  trajectory  was  small.  If  the  resolution  of  the  ocean  velocity  is  high  in 
one  grid  cell  of  the  Regional  NCOM  data,  the  trajectory  could  be  estimated  more 
precisely.  Some  NCOM  prediction  data  sets  have  higher  resolution,  less  than  one 
kilometer.  Thus,  research  would  be  able  to  attempt  to  detennine  whether  an  optimal 
trajectory  can  be  estimated  by  using  both  higher  resolution  data  and  our  correction 
methods.  If  the  two  conection  methods  can  be  applied  to  this  situation,  operators  would 
get  better  trajectory  estimation. 

Finally,  this  study  only  considered  the  ocean  cunent  effects.  The  position  of  an 
underwater  glider  is,  however,  determined  by  both  vehicle  and  environment  models  as  we 
discussed  in  Chapter  II.  Thus,  it  is  expected  that  a  more  precise  trajectory  can  be 
estimated  by  considering  both  the  ocean  current  and  a  glider’s  vehicle  model  effects. 
Some  previous  research  studied  the  glider’s  dynamic  model  effects  without  considering 
the  ocean  current  effects,  but  more  expansive  research  would  incorporate  the  two  effects 
simultaneously.  Even  though  this  recommended  future  research  is  complex  and  requires 
more  time,  it  would  be  able  to  estimate  a  more  precise  trajectory. 
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APPENDIX  A:  RESULTS  OF  TRAJECTORY  ESTIMATION 


A.  PHASE  ONE:  DR  AND  DRNCOm  TRAJECTORIES 
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Phase  One:  Dive  #10 
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Phase  One:  Dive  #12 
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Phase  One:  Dive  #15 
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Phase  One:  Dive  #19 
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Phase  One:  Dive  #25 
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Phase  One:  Dive  #34 
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Phase  One:  Dive  #37 
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Phase  One:  Dive  #43 


Phase  One:  Dive  #44 


Phase  One:  Dive  #45 
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Phase  One:  Dive  #46 


Phase  One:  Dive  #47 


Longitude 


Phase  One:  Dive  #48 
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Phase  One:  Dive  #49 
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Phase  One:  Dive  #52 
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Phase  One:  Dive  #55 
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Phase  One:  Dive  #61 
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Phase  One:  Dive  #63  Phase  One:  Dive  #64 
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Phase  One:  Dive  #65 
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Phase  One:  Dive  #68 
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B.  STEP  ONE:  DR  TRAJECTORIES  WITH  DEPTH-AVERAGED 
CORRECTION  METHOD 


Step  One:  Dive  #1 


Step  One:  Dive  #2 


Longitude 


Longitude 
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Step  One:  Dive  # 7 


Step  One:  Dive  #8 


Step  One:  Dive  #9  Step  One:  Dive  #10 
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Step  One:  Dive  #12 
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Step  One:  Dive  #13 


Step  One:  Dive  #14 
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Step  One:  Dive  #18 
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Step  One:  Dive  #19 


Step  One:  Dive  #20 
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Step  One:  Dive  #25 


Step  One:  Dive  #26 
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Step  One:  Dive  #30 


91 


Latitude  Latitude  Latitude 


Step  One:  Dive  #31 


Step  One:  Dive  #32 


Longitude 


Step  One:  Dive  #33 


Step  One:  Dive  #34 


Step  One:  Dive  #35 


Step  One:  Dive  #36 


92 


Latitude  Latitude 


Step  One:  Dive  #37 
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Step  One:  Dive  #43 
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Step  One:  Dive  #49 
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Step  One:  Dive  #54 
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Step  One:  Dive  #55 
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Step  One:  Dive  #61 
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Step  One:  Dive  #66 
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Step  One:  Dive  #67  Step  One:  Dive  #68 
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STEP  TWO:  DR  TRAJECTORIES  WITH  DEPTH-DEPENDENT 
CORRECTION  METHOD 


Step  Two:  Dive  #1 


Longitude 


Step  Two:  Dive  M2 


Longitude 


Stop  Two:  Dive  MS 


Stop  Two:  Dive  #6 


Longitude 


99 


Latitude  Latitude  Latitude 


Step  Two:  Dive  #7 
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Step  Two:  Dive  #13 


Step  Two:  Dive  #14 


Stop  Two:  Dive  #15 
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Step  Two:  Dive  #19  Step  Two:  Dive  #20 
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Step  Two:  Dive  #23  Step  Two:  Dive  #24 
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Step  Two:  Dive  #25  Step  Two:  Dive  #26 
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Step  Two:  Dive  #31 
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Step  Two:  Dive  #37 
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Step  Two:  Dive  #43 
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Step  Two:  Dive  #49 
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Step  Two:  Dive  #55 
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Step  Two:  Dive  #61 
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Step  Two:  Dive  #67 
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STEP  THREE:  ESTIMATED  OPTIMAL  TRAJECTORIES  FROM 
C2  BY  USING  ITERATION  METHOD 


Step  Three:  Dive  #1 
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Step  Three:  Dive  #7 
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Step  Three:  Dive  #13 
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Step  Three:  Dive  #25 
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Step  Three:  Dive  #31 
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Step  Three:  Dive  #37 
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Step  Three:  Dive  #38 


Step  Throe:  Dive  #39 
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Step  Throe:  Dive  #40 


Step  Three:  Dive  #41  Step  Three:  Dive  #42 
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Latitude  Latitude  Latitude 


Step  Three:  Dive  #43 


Step  Three:  Dive  #44 


Step  Thrco:  Dive  #45 
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Step  Throe:  Dive  #46 


Step  Three:  Dive  #47 


Step  Three:  Dive  #48 
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Latitude  Latitude  Latitude 


Step  Three:  Dive  #49 


Step  Three:  Dive  #50 
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Step  Three:  Dive  #51 


Step  Three:  Dive  #52 


Step  Throe:  Dive  #53 


Step  Three:  Dive  #54 
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Latitude  Latitude  Latitude 


Step  Three:  Divo  #55 


Step  Three:  Dive  #56 
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Step  Three:  Divo  #57 
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Latitude  Latitude  Latitude 


Step  Three:  Dive  #61 


Step  Three:  Dive  #62 
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Step  Three:  Dive  #63 


Step  Three:  Dive  #64 


36  6836 


36  6832 


36.6828 


36  6824 


-121  9344  -121  934  -121  9336 

Longitude 


36  682 

36.6816 

o 

•a 

3 

8 

to 

j 

36  6812 
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Stop  Three:  Divo  #65 


Step  Three:  Dive  #66 


Longitude 


121 


Latitude 


Step  Three:  Dive  #67 


Step  Three:  Dive  #68 
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Latitude  Latitude  Latitude 


E.  STEP  FOUR:  COMPARISON  THREE  CORRECTED 

TRAJECTORIES 


-121  931  -121  9308 
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Step  Four:  Dive  #3 
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Step  Four:  Dive  #2 
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Step  Four:  Dive  #5 


Step  Four:  Dive  #6 
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Latitude  Latitude  Latitude 


Step  Four:  Dive  #7 


Step  Four:  Dive  #8 


Step  Four:  Dive  #9 


Stop  Four:  Dive  #10 
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Latitude  Latitude  Latitude 


Step  Four:  Dive  #13 


Stop  Four:  Dive  #14 
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Stop  Four:  Dive  #15 


Stop  Four:  Dive  #16 
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Stop  Four:  Dive  #17  Step  Four:  Dive  #18 


125 


Latitude  Latitude  Latitude 


Stop  Four:  Dive  #19 


Stop  Four:  Dive  #20 


Stop  Four:  Dive  #21 


Step  Four:  Dive  #22 


Stop  Four:  Dive  #23  Stop  Four.  Dive  #24 
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Latitude  Latitude  Latitude 


Stop  Four:  Dive  #25 


Stop  Four:  Dive  #26 


Stop  Four:  Dive  #27  Stop  Four:  Dive  #28 


Step  Four:  Dive  #29 


Step  Four:  Dive  #30 
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Latitude  Latitude  Latitude 


Step  Four:  Dive  #31 


Step  Four:  Dive  #32 


Longitude  Longitude 


Stop  Four:  Dive  #33 


Step  Four:  Dive  #34 


Stop  Four:  Dive  #35  Step  Four:  Dive  #36 
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Latitude  Latitude  Latitude 


Step  Four:  Dive  #37 


Longitude 


Step  Four:  Dive  #38 


Step  Four:  Dive  #39  Step  Four:  Dive  #40 


Step  Four:  Dive  #41  Stop  Four:  Dive  #42 
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Latitude  Latitude  Latitude 


Stop  Four:  Dive  #43 


Step  Four:  Dive  #44 
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-121  9392  -121  9391  -121  939 
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Step  Four:  Dive  #47 
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Latitude  Latitude  Latitude 


Stop  Four:  Dive  #49 


Stop  Four:  Dive  #50 
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Latitude  Latitude  Latitude 


Stop  Four:  Dive  #55 


Step  Four:  Dive  #56 
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Stop  Four:  Dive  #61 


Stop  Four:  Dive  #62 
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-121  9352  -121.9348  -121  9344  -121.934 

Longitude 


Stop  Four:  Dive  #63 


Longitude 
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Step  Four  Dive  #65  Stop  Four:  Dive  #66 


133 


Latitude 


Step  Four:  Dive  #67 


Step  Four:  Dive  #68 


Longitude 
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APPENDIX  B:  COINCIDENCE  OF  DR  AND  DRNCOm  TRAJECTORY 
WITH  THE  DEPTH-DEPENDENT  CORRECTION  METHOD 

In  this  appendix,  we  prove  the  results  coincidence  of  DR  (Ti)  and  DRNCom  (T2) 
trajectory  with  the  depth-dependent  correction  method  (see  Figure  27  and  Figure  28, 
respectively).  It  helps  us  save  some  effort  and  process  efficiently  for  improving  our 
method  in  this  study.  We  form  a  hypothesis  that  both  corrected  trajectories  with  the 
depth-dependent  correction  method  have  a  same  result.  By  proving  this  hypothesis,  the 
only  trajectory  that  needs  to  be  corrected  would  be  the  Ti  to  develop  our  method. 


(X*  yf) 


Figure  27.  Lateral  View  of  the  Eastward  Component  of  the  T i  with  the  Depth- 

dependent  Correction  Method. 


Let  there  be  N  steps  from  start  to  end  of  each  trajectory.  Then,  the  final  position 
of  T i  and  T3  can  be  expressed  as  P dr  (  Xdr ,  Ydr )  and  Pgps  ( X N ,  YN  ),  and  the  distance  error 
between  Pdr  and  Pgps  also  can  be  presented  as  AXW  =  XN  -  Xdf ,  AY"'  =YN—  Ydr .  Now, 

let  us  correct  the  Ti  on  the  x-direction  by  applying  the  depth-dependent  correction 
method  as  shown  in  Chapter  III. 
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The  corrected  water  displacement  ( Ax" )  at  each  time  ( tn )  can  be  expressed  as 
shown  in  Equation  (34)  by  applying  the  Equation  (27)  and  Equation  (32)  in  Chapter  III. 


y  fz  (‘At- 

/  jJ  i  Ji  i  . ncom  /  ~,ncom((\  +\ 

Ax:  =  AX~^ - .where  fr  =  <Z"A  />"■ 

YJ*  ft  ^  Ui  vV>  ^i)  Wmax 

i=l 

and  wmax  =  max  (0,t/)],  /  =  1,2,-,  N 


(34) 


,  ncom 
N 


Figure  28.  Lateral  View  of  the  Eastward  Component  of  the  T?  with  the  Depth- 

dependent  Correction  Method. 


At  this  time,  we  correct  the  T2  on  the  x-direction  by  applying  the  depth-dependent 
correction  method.  Let  the  final  position  of  T2  be  Pncom  ( X'(:om ,  Y(com  ),  and  the  distance 

error  between  Pncom  and  Pgps  also  be  AXnw  =  XN  -  X(om ,  A  Y'lw  =  YN-  Y(com  .  Moreover, 
with  the  relationship  of  position  between  Ti  and  T2,  the  position  of  T2  at  each  time  ( tn ) 


can  be  expressed  as  follows. 


\rncom  \rdr  ,  \  1  ncom  /  -  \ 

Xn  =Xn  +2^Ui  (ZP  tPAti 
1=1 

Y™  =Yndr  +fjvrm(zi,ti)Ati 


i= 1 


(35) 
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Then,  Equation  (35)  can  be  presented  as  shown  in  Equation  (36). 


n  n  ,,ncom(7  t  \  1inCOm(()  t  ) 

\rncom  _  \rdr  .  \  '  ncom  /  +\\+  _  vdr  .  \  1  a, 

+2-fM/  yZi^i)Ati  ~An  +UmaxZj^~ - A^ 


i=l 


ncom  /r\  ,  \ 

i= 1  W(  (<U,)  Mmax 


=  where  =  max  [„”  (0,  <,)],<  =  1,2,--,JV 


(36) 


In  addition,  the  final  position  can  be  presented  as  shown  in  Equation  (37). 

xr=K+“^tfif:  Ati 

i= 1 

Consequently,  the  total  corrected  water  displacement  will  be 


(37) 


AX”  =  XN-  X~  =  AX"  X 


(38) 


i=l 


The  corrected  water  displacement  ( Ax"w )  at  each  time  ( tn  )  can  be  expressed 
similar  to  the  Equation  (34)  by  substituting  Equation  (38)  into  Equation  (34)  as  shown  in 
Equation  (39). 


A<"  =  AX”  A! - =  (AX' -f - 


a<, 


1=1 


i=l 


AX" 


—  u 


i= 1 

=  Axw-u  y  fzflAt 

n  max  J  i  J  i  \ 


(39) 


i= 1 


Finally,  the  corrected  position  of  C2  and  C5  (^xfr  ,Xcn',com)  at  each  time  is  the 
same,  and  it  can  be  expressed  as  shown  in  Equation  (40). 


x:dr  =xdj+Ax: 


x ; 


cncom  _  ncom  \\nw  _ 


A  f 


Xdr+u  Yfzf‘At.  +  Axw-u  Y  rf‘At. 

n  max  J 1  J 1  1  n  max  J 1  J 1  1 


i= 1 


j  v 


i= 1 


(40) 


=xf+Axr =x 


cc/r 
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APPENDIX  C:  ITERATION  METHOD  FOR  ESTIMATING  AN 
OPTIMAL  UNDERWATER  TRAJECTORY  OF  GLIDER 


This  appendix  describes  the  iteration  method  to  estimate  an  optimal  trajectory  by 
applying  the  depth-dependent  correction  method  with  C2.  C2  is  the  corrected  trajectory 
using  the  depth-dependent  correction  method  with  DR  trajectory.  In  this  study,  the  real 
underwater  trajectory  of  glider  is  unknown,  thus  an  optimal  trajectory  development  is 
needed  to  determine  which  corrected  trajectory  is  more  precise  by  comparing  the  two. 
The  iteration  method  is  described  mathematically  as  follows. 

Let  R(Y)  be  the  glider’s  position  vector  with  the  function  of  time,  U  (R,*)  be  the 
glider’s  velocity  vector  relative  to  the  ocean  currents  with  the  function  of  the  glider’s 
position  and  time,  and  Uncom(R ,i)  be  the  water  velocity  vector  with  the  function  of  the 
glider’s  position  and  time. 

Generally,  the  glider’s  position  at  time  tn  can  be  expressed  as  shown  in  Equation 

(41). 

R(0  =  R(f0)  +  f"  U(R,  t)dt  =  R(f0)  +  f"  Ug(R,0*  +  f"  Uw(R,0  dt  (41) 

J  t0  J  tQ  J  tQ 

In  this  study,  the  glider’s  DR  position  does  not  consider  the  water  velocity,  so  its 
DR  position  would  be 

R*(0  =  R(*o)+ f”  U  (R,0<*  (42) 

J  t0 

Then,  the  final  position  of  glider  at  time  tN  can  be  presented  as  shown  in  Equation 
(43)  by  substituting  Equation  (42)  into  Equation  (41). 

R«»)  =  R*«»)+l'"U™(R.O^  (43) 

The  C2  is  corrected  from  the  Ti  by  applying  the  depth-dependent  correction 
method  against  the  distance  difference  between  Pdr  and  Pgps  (AX"' ,  see  Figure  4).  Now, 
we  develop  the  x-dircction  component  of  the  C2.  Let  the  glider’s  x-dircction  component 
of  final  position  vector  at  time  tN  be  X  (tN). 
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Then,  it  can  be  expressed  as  follows. 

•'*0 

where  (44) 

Xdr(tN)  :  x -direction  component  of  DR  position  vector 

uncom(R,t) :  Eastward  component  of  the  Regional  NCOM  velocity  vectors 

Then,  the  distance  difference  between  Pncom  and  Pgps  ( AX""' )  would  be 

AX'  =  =  f"  u„(K,t)dt  (45) 

•'*0 

To  correct  AX"'  ,  Equation  (46)  is  obtained  by  applying  Equation  (29)  and 
Equation  (30). 

N 

AX"  =  C£f;f‘At, 

i= 1 

(46) 

C  =  — - ,  C:  Correction  Coefficient 

tf'fAi, 


In  addition,  the  distance  difference  between  T2  and  T3  (Ax™)  at  each  time  ( tn ) 
would  be 

Y^f-flAt, 

Ai-;  =  At,  =  AX~  (47) 

Tfif'  A/, 

1=1 

After  substituting  Equation  (47)  into  Equation  (40)  in  Appendix  B,  the  eastward 
component  of  glider’s  corrected  position  vector  Xcdr(tn)  at  time  tn  would  be 


Zf'fAt, 

X“"(t„)  =  XJr(t,)+ AZ"4=! -  (48) 

Tf’fAt, 


However,  since  the  variables  in  Equation  (48),  fz ,  and  f‘  depend  on  the  glider’s 
position  and  time,  the  equation  is  non-linear.  Thus,  it  is  difficult  to  get  an  exact  solution. 
Therefore,  we  use  the  iteration  method  to  reach  close  to  the  exact  solution. 
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Let  the  glider’s  eastward  component  of  position  vector  after  iteration  at  time  ( tn ) 
be  Xm+1(tn),  then 

tfrimtfim)  At, 

Xm+1  (tn )  =  X*  )  +  AXm'  -M - ,  m  =  1,  2, . . . ,  M  - 1,  M 

;=i 

where  f‘(m) :  Factors  on  Xm,  Ym 

To  determine  how  many  times  iteration  will  be  done,  we  need  to  set  a  criterion. 
We  set  the  relative  error  of  10-5  for  this  study.  Let  the  reference  and  iterated  trajectory  be 
Tm  and  Tm+i,  the  sum  of  distance  difference  between  two  trajectories  be  Dm  ,  and  the  total 
length  of  the  reference  trajectory  be  Lm  at  each  step.  Then  the  relative  error  ( Erd  )  at  each 
step  can  be  calculated  by  dividing  Dm  by  Lm  .  Equation  (50)  presents  this  procedure 
mathematically. 

Erel=Dm/Lm,rn  =  l,2,-,M-l,M 

N 

A,=Zor‘-*.') 

n= 1 

where  (50) 

Xm+l(tn )  :  The  eastward  position  of  glider  on  the  iterated  trajectory  at  time  tn 
Xm(tn) :  The  eastward  position  of  glider  on  the  reference  trajectory  at  time  tn 

If  Erel  is  less  than  10-5,  the  iteration  is  allowed  to  be  tenninated. 
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